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Abstract 

We are entering an era where a number of large-scale lattice simulations of 
four-dimensional supersymmetric theories are under way. Moreover, proposals 
for how to approach such studies continue to progress. One particular line 
of research in this direction is described here. General actions for super-QCD, 
including counterterms required on the lattice, are given. We obtain the number 
of fine-tunings that is required, once gauge and flavor symmetries are accounted 
for, provided Ginsparg- Wilson fermions are used for the gauginos. We also 
review and extend our recent work on lattice formulations of N = 4 super- Yang- 
Mills and M = 1 super- Yang-Mills that exploit Ginsparg- Wilson fermions. 
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1 Introduction 



1.1 Motivations 

Supersymmetric lattice field theories are motivated by the desire to obtain nonper- 
turbative information that cannot be obtained by other means. Having reviewed the 
aims of lattice supersymmetry recently in [1], here we will offer only a brief reiteration, 
adding a few remarks having to do with recent developments. 

Definition. A classic example of a nonperturbative anomaly is the Witten 
anomaly [2]. A more recent example is [3]. It has been argued that a nonperturbative 
supersymmetry anomaly exists [4]. The lattice approach is a tool to investigate this 
question. 

Nonholomorphic quantities. In a theory of chiral superfields <j), holomorphic 
quantities w((f>) are protected by nonrenormalization theorems. When combined with 
symmetries and "the power of holomorphy" [5], much can be learned about strong 
dynamics in supersymmetric theories. Not so for nonholomorphic quantities k((f>,4>); 
little is known about them in the strongly coupled regime. The supersymmetry- 
breaking soft-terms that determine spectra and couplings in supersymmetric exten- 
sions to the Standard Model depend on the nonholomorphic Kahler potential; see 
for example [6,7]. In fact, it has recently been realized that strong hidden sector 
effects can lead to significant modifications of the observable sector soft terms [8], 
with the potential to solve some long-standing phenomenological problems, such as 
the [i/Bp problem in gauge mediation models [9]. A strong coupling computational 
method is needed in order to say anything definite. To answer the crucial question 
in [9] — the sign of an anomalous dimension — lattice artifacts at the level of 10-20% 
may be tolerable. One goal of the research that we report here, in £J2} is to develop 
lattice super-QCD as a tool to study this sort of problem. 

Dynamical supersymmetry breaking. New strong gauge theory interac- 
tions are commonly employed to split the superpartners from the observed Standard 
Model spectrum^ Strongly coupled messenger sectors and compositeness within the 
supersymmetric standard model [13-18] can provide both economy to the models and 
difficult strong coupling questions of phenomenological importance at scales of a few 
TeV, hence relevant to the Large Hadron Collider [19]. Moreover, such models are 

Reviews include [10,11] for supergravity mediation of gaugino condensation, and [12] for gauge 
mediation. 
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well- motivated by warped string compactifications such as have been explored in [20], 
or the older, perturbative string compactifications on toroidal orbifolds [21-24] whose 
phenomenology has been extensively studied [25-32]. Advances in lattice supersym- 
metry move us toward addressing these questions. Indeed, it is important to continue 
exploring supersymmetric models other than the most popular scenarios such as the 
constrained minimal supersymmetric standard model (CMSSM), although new possi- 
bilities in these well-studied models continue to be uncovered, as in [33] . Some of the 
most promising models of dynamical supersymmetry-breaking involve chiral gauge 
theories. This is a very difficult problem that we will not address here, though there 
is some interesting recent work on formulations using Ginsparg- Wilson fermions and 
strong Yukawa couplings [34-38] . 

Gauge-gravity duality. Recently lattice methods have contributed [39-41] to 
the evolving understanding of the relationship between supersymmetric gauge theory 
and theories of quantum gravity; in particular, string/M-theory. In fact, the a' and 
string loop corrections to certain effective supergravity descriptions of string theory 
in nontrivial backgrounds are supposed to be encoded in corrections to the 't Hooft 
limit of the gauge theory. Obviously, the lattice theory at finite iV and coupling g will 
capture these effects, though we must still take the continuum limit — which includes 
somehow restoring the supersymmetry broken by the lattice regulator. 

It is worth mentioning that in the case of matrix supersymmetric quantum me- 
chanics, a non-lattice approach has been developed with considerable success in [42] 
and subsequent articles. These authors fix the gauge and work directly in momentum 
space with a sharp cutoff, for this case of 1 + dimensions. Supersymmetry breaking 
by the regulator is believed to be especially mild. Certainly it vanishes as the mo- 
mentum cutoff is removed to infinity, since one then obtains the unregulated theory, 
which is finite and requires no subtractions. 

A background independent, nonperturbative formulation of superstring theory is 
not known. Nevertheless, the theory is in much better shape due to successes that 
address nonpertubative and background dependent aspects: M(atrix) theory [43,44], 
the AdS/CFT correspondence [45-47], the PP-wave limit [48-51], etc. It would be 
very interesting to study these formulations through their relation to super- Yang- 
Mills (SYM). The Matrix theory formulations of string/M-theory, and the AdS/CFT 
correspondence, are expressed in terms of quantum theories of dimensionally reduced 
SYM. The vacuum of the gauge theory is believed to have a gravitational mean- 
ing. Detailed studies of the SYM vacuum might provide useful information with a 
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gravitational interpretation. 
1.2 Challenges 

Dondi and Nicolai [52] pointed out long ago that since the supersymmetry algebra 
closes on the generator of infinitesmal spacetime translation, which is explicitly broken 
by the discretization, the supersymmetry algebra invariably must be modified on the 
lattice. The obvious option is to have it close on discrete translations. Since the 
Leibnitz rule does not hold on the lattice, the supersymmetry algebra will be violated 
for general polynomials of lattice fields; interacting supersymmetric theories will not 
be invariant with respect to the lattice supersymmetry [53,54]. 

The non-invariance of the interacting lattice action is an 0(a) effect (a is the 
lattice spacing) that disappears if one takes the continuum limit of the lattice ac- 
tion. Unfortunately, in the quantum theory these violations, which correspond to 
0(a) irrelevant operators that supplement the continuum action, play off against ul- 
traviolet (UV) divergences to give infinite violations of supersymmetry in the a — > 
limit. Another way of stating it is this: non-irrelevant supersymmetry- violating op- 
erators allowed by the symmetries of the lattice action will be radiatively generated 
(e.g., a mass term for scalar partners of the gauge boson in extended SYM theo- 
ries); supersymmetry- violating relative renormalizations of terms already present in 
the bare action also will occur (e.g., quartic scalar self-couplings with a coefficient 
other than g 2 /2, where g is the gauge coupling). 

In fact, the situation is similar to the chiral limit for Wilson fermions, where the 
bare mass must be fine-tuned in order to cancel the effects of the 0(a) suppressed 
irrelevant Wilson mass operator. Fine-tuning of counterterms to achieve the de- 
sired continuum limit is in principle always possible, but an efficient nonperturbative 
method is required if one wishes to address strongly interacting theories. This is the 
sort of approach advocated in $2]-£|I]of this review. 

A contrasting situation is the one where symmetries of the lattice theory over- 
come the difficulty of supersymmetry-violating renormalizations. An example is the 
domain wall fermion formulation of M = 1 SYM discussed in $51 In that case, lattice 
chiral symmetry in the form of Ginsparg- Wilson (GW) fermions [55] prevents addi- 
tive renormalization of the gluino mass in the continuum limit, and hence the only 
supersymmetry-violating non- irrelevant operator is forbidden in that limit by setting 
the bare mass to zero [56-64]. 

In fact, use of symmetries to prevent bad renormalization is a much more general 



3 



approach that can be applied in a number of theories [65, 66] including some that 
are four- dimensional [67-69]. Some early detailed studies of these types theories have 
appeared recently [70-75]. Many aspects of this approach were discussed in recent 
reviews [1,76,77], which contain more extensive references to this line of research. We 
will not dwell on such formulations. 

Since M = 4 SYM and super-QCD have scalars, even with lattice chiral symmetry 
for the fermions the scalar masses and couplings (Yukawa and quartic) will receive 
divergent non-supersymmetric corrections in the continuum limit and must be (non- 
perturbatively) fine-tuned to supersymmetric values. For this reason these theories 
have always seemed impractical by the fine-tuning approach! 

We have argued in [78] that for M = 4 SYM this opinion is overly pessimistic. In 
this review we extend our line of reasoning to super-QCD (SQCD) where the situation 
is, confessedly, more challenging Jf| First, however, let us concentrate on the reasons 
why M = 4 is more practical that might be naively concluded: 

1. If one uses GW fermions the gluinos can be kept massless. 

2. The SU(4)r symmetry can be preserved, restricting renormalizations. 

3. The parameters that must be tuned consist of: 

• one scalar mass, 

• two/four quartic couplings (N c = 2,3 vs. iV c > 3), and 

• one Yukawa coupling. 

Next we summarize aspects of the fine-tuning procedure that help to make nonper- 
turbative adjustment of four/six parameters "practical." 

The Yukawa coupling can be tuned by rescaling the scalar kinetic term. This 
is obvious because the Yukawa coupling strength y can always be absorbed into a 
redefinition of the scalar fields <fi —>■ <f)/y, causing it to reappear in the scalar kinetic 
term. Thus, all tunings can be done by adjusting bosonic terms in the action. This 
allows the tunings to be done by the "Ferrenberg-Swendsen method" [81-83] , explor- 
ing a wide swath of coupling constant space "offline" from the results of a single 
Monte-Carlo simulation. The parameter range available with good statistics can be 

2 Other approaches to TV" = 4 SYM include Refs. [67-69], which involve orbifold or twisted super- 
symmetry lattices. 

3 Other approaches to SQCD have been explored in [79,80], albeit in a two-dimensional context. 
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enlarged using multicanonical techniques [84-89]. Thus we arrive at the encourag- 
ing result that all fine-tuning can be performed through an offline analysis; i.e., new 
simulations that require large numbers of fermion matrix inversions during molecular 
dynamics trajectories are not required. 

The message is this: one need only generate a set of configurations that coarsely 
cover the parameter space in the vicinity of the fine-tuned lattice parameters. In 
practice, this neighborhood of the Af = 4 SYM or SQCD point in parameter space 
would be determined by starting on very small lattices, and in fact perturbative 
calculations would give us a good idea where to begin for weak bare couplings on 
such lattices. This is because on a small lattice there is not much separation between 
the UV and infrared (IR), and hence the effective coupling remains weak at the IR 
scale. Simulations can then be used to move into stronger coupling regimes, so that 
one bootstraps upon previous results in order to stay in the supersymmetric window 
for bare couplings. Modest computational resources would be able to perform all the 
offline fine-tunings, and to carefully study the location of the supersymmetric point 
in the bare lattice parameter space. 

Similar statements hold for Af = 2 SYM theory, though we do not go into details 
here. However, we note that fine-tuning with Wilson fermions has been analyzed by 
Montvay in [90]. By contrast, if GW symmetry is exploited, the SU(2)r symme- 
try of the continuum can be preserved, which reduces the number of counterterms 
significantly. 

1.3 Overview 

We begin in $2] with some new results, providing the most general continuum La- 
grangians for SQCD consistent with the symmetries that the lattice will preserve. 
This then allows us to write down the lattice actions including all counterterms that 
must be included in order to fine-tune to the supersymmetric point. We enumerate 
the number of fine-tunings in each case, and show that they can be accomplished 
through entirely bosonic reweighting, with one exception. In £J3] we review the re- 
sults of our previous study of similar approach in Af = 4 SYM, extending some of 
the discussion of how tuning can be implemented. In §|4] we give a discussion of the 
multicanonical reweighting method, and explore how it can be implemented for the 
theory described in §|3j In £J5] we describe recent large-scale simulations of Af = 1 
SYM with a domain wall fermion implementation, adding some results on nonlinear 
chiral extrapolations. The domain wall fermion implementation of Af = 1 SYM is 
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in the spirit of §|2|-£j3j except that the GW symmetry is in this case powerful enough 
to prevent any counterterms that would have to be fine-tuned. Conclusions and an 
Appendix follow these main sections. 



2 Lattice super-QCD 

Super-QCD extends QCD by adding a fermionic partner for the gluon and scalar 
partners for the quarks. Since the theories that we are interested in here are new 
gauge interactions that are strong at scales of a TeV or greater, it is only related to 
QCD by way of analogy. The gauge groups that we consider here will be SU(N), 
and the number of flavors will be Nf. The continuum theory is briefly reviewed in 
Appendix [A] 

2.1 SU(2) gauge theory 

The gauge sector consists of the gauge boson and the gaugino A. The gauge action 
will be formulated using a massless GW fermion A and the Wilson plaquette action or 
some improved version of it. For numerical stability, it will be necessary to simulate 
at nonzero gaugino mass and then extrapolate to chiral limit ni\ — > 0. 

The matter sector of the theory consists of Pj , Qi, 1 = 1,..., Nf chiral superfields, 
each containing a complex scalar and a left-handed Weyl fermion. Here we distinguish 
them based on U(l)v charge, +1 for P, —1 for Q, a symmetry of the continuum theory 
that will be preserved exactly in the lattice formulation. On the other hand, Pj and 
Qi are both fundamentals of SU(2), transforming identically. We denote scalar and 
left-handed Weyl fermion components by pi, qi and Xpi-> Xqi respectively. 

2.1.1 General form of invariants 

The scalar quadratic SU(2) invariants are given in Table [TJ We will impose two other 
constraints that are symmetries of the continuum Lagrangian, given in the Appendix, 
Eq. IA.1I The first is CP conservation, and hence reality of the coefficients. The second 
is a Z 2 exchange symmetry that we will call S: 



Then the most general mass term for the scalars, suppressing flavor indices, is: 



S : Vi *-> q T , p} «-> q*!, X P i «-> Xqi, X P i <-> Xqi- 



(2.1) 



m 



2 (p'p + q^q) ■ 



(2.2) 
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SU(2) invariant U(l)y charge 


SU(2) invariant U(l)y charge 


ptp 
p T eq 


gig 
p+eg* 


q^p 2 
p T ep 2 
qhq* 2 


p^g -2 
p^ep* -2 
g T eg -2 



Table 1: Bilinear SU(2) invariants, flavor indices suppressed. 



The other SU(2) x U(l)y invariant mass ter 

m' 2 (p T eq + p^eq*) (2.3) 

is ruled out by the S exchange symmetry (12.11) . Later, for number of flavors Nf > 1, 
we will impose SU(Nf) p x SU(Nf) q flavor symmetry constraints, and specify the 
corresponding matrix structure of m 2 . This flavor symmetry also forbids (12.3H . 

The independent quartic invariants, only taking into account SU(2) gauge invari- 
ance and U(l)y at this point, are built from the bilinear invariants in the Table [U 
Combining the U(l)y neutral bilinears, and imposing CP and S symmetries, we 
arrive at the "(0,0)" quartic Lagrangian: 

Ao,o) = HiP ] P? + (<? f 9) 2 ] + A 2 pW<7 + A 3 |peg| 2 + K[{ptqf + (pV) 2 ]- (2.4) 

The other (0,0) term that can be obtained from the quadratics in the table is 

\ h (p ] p + q ] q)(ptq + p*eq*) ) (2.5) 

but it is ruled out by the S symmetry (12.11) . Likewise we combine the U(l)y charged 
bilinears to obtain the "(2,-2)" quartic Lagrangian: 

£(2,-2) = v x \q ] p\ 2 + v 2 (\p T tp\ 2 + \q T tq\ 2 ) 

+z/ 3 [(p T ep)(g T eg) + (p t ep*)(g t eg*)] 

+z/ 4 [(g t p)(g T eg + p f ep*) + (p t g)(g t eg* +p T ep)]. (2.6) 

In fact, the v± term will violate the nonabelian flavor symmetry SU(Nf) p x SU(Nf) q 
and we will end up discarding it for all but the one flavor case. 

4 Here and below T denotes transpose, while e is the two-dimensional Levi-Cevita tensor with 
convention ei2 = 1 = —621. 
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Finally, note that we have eliminated other SU(2) invariant quartic operators 
through relations such as 



(p\v a pj)(p> a pL) 



2 {pW*k){pj£Pl) 



(2.7) 



using the Fiertz identity 
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05 - - 



(2.8) 



0=1 



In Eq. (12. 7ft . I,J,K,L are flavor indices. Identities like ( 12. 7ft will relate the general 
Lagrangian that we are writing down to the supersymmetric theory, since in the 
latter the quartic interactions are typically expressed in the form of the left-hand side. 
We will return to this below, once flavor symmetry constraints have been taken into 
account. For now we merely state that the D-term Lagrangian in the supersymmetric 
theory is: 



On comparing (12. 9p to the quartic interactions of the general theory, (12.41) and 
(12. 6p . we see that we have many more quartic interaction parameters than in the 
continuum theory, where there is only one type of term with a strength determined 
by the gauge coupling. It will be seen below that this is a general feature of the SQCD 
theories: many fine-tunings are needed due to a large number of quartic couplings 
that are allowed. We postpone the precise count of finely-tuned parameters until we 
take into account flavor structures. We will do that shortly, but first we complete our 
general parameterization by considering the Yukawa couplings. 

Here it is not hard to check that the SU(2) gauginos A a , a = 1,2,3 give rise to 
the following unique CP and S symmetric, hermitian Yukawa interactions, written in 
two- component notation: 




(2.9) 




yilp^Xp + Xp^p + q^Xg + x T q M] 

+y 2 [p T e\x q - Xq^P + q T ^X P - xjAeg] 



(2.10) 



where 



A = X a t' 



A = X a t' 



t a = -a' 
2 



.a 



(2.11) 
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and a a are Pauli matrices. It is obvious that this is hermitian. CP acts on the fields 
according to: 



P > P ) Xp > Xpi Q *9 i Xq > Xq: 

X _> A* = X a t aT = X T , X — > A T . (2.12) 

It can be checked that ( 12.101) is invariant under this symmetry. In doing this one must 
keep in mind the rules of two-component spinors as it relates to Grassmann fields, 
such as 

X a x P = Kxt = -xtK = X P *\ a& = X P A a , (2.13) 

see Appendix B of [91]. In the supersymmetric target theory, y 2 = and y± = V2g, 
where g is the gauge coupling, together with the field redefinition A — ► iX, X — > — iX. 
See Appendix $A] for further details. 

The 1/2 terms are forbidden if there is more than one flavor, which is quite useful 
since directly tuning fermionic interaction terms is most likely not practical. Two 
approaches to the tuning of the y\ term will be discussed below, both of which involve 
tuning bosonic terms relative to this fermionic term. 

2.1.2 One flavor 

The flavor symmetry is simple to analyze when Nf = 1, since there are no flavor 
indices to add to the operators that we have just written down. All but the first of 
the (2,-2) operators in (12. 6p vanish identically, eliminating ^2,3,4 from consideration. 
All of the (0,0) operators in (12 ,4p are allowed, and the mass term (12.21) and both types 
of Yukawa terms ( 12. 10}) are also unrestricted. Altogether nine parameters must be 
fine-tuned: 

• one scalar mass m 2 , 

• six quartic couplings Ai^.3,4,5, v\, and 

• two Yukawa couplings yip- 

As mentioned already, in the Nf = 1 theory two "fermionic" parameters that must 
be fine-tuned, y\^. In a simulation they are buried in the fermionic determinant, as 
far as the Boltzmann weight in configuration space is concerned. Unlike the bosonic 
parameters, they cannot be adjusted "offline" by reweighting techniques. If we appeal 
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to the method advocated in [78], then one of the Yukawa fine-tunings can be done 
equivalently through introducing a p, q scalar field strength Z<p in front of the kinetic 
term: 

Z^\Dp\ 2 + \Dq\ 2 ). (2.14) 

This rescaling of the scalars at tree level then provides a lever to adjust y\. (Of course, 
the mass must also be rescaled m 2 — > m' 2 pa Z^m 2 , and then more precisely tuned, 
to keep the theory near the desired physical mass.) What is needed, therefore, is a 
lattice symmetry that enforces yi = 0. Unfortunately, none seems to be available. 
The chiral symmetry 

p- e *p, q^e l9 q, X P ^e ld Xp , Xq ^e w Xq (2.15) 

is anomalous and is of no help. We conclude that the Nf = 1 SU(2) SQCD will 
be very difficult to study in the current formulation, due to the additional type of 
Yukawa coupling. 



2.1.3 Two flavors 



The target theory has a SU(2) p x SU(2) q flavor symmetry that we will preserve in 
the lattice action. The unique mass term arising from (12. 3p is 



m 



(2.16) 



7=1 



with Nf = 2. The quartic terms in (12.41) allow for six flavor symmetric terms 

aS 1} (p\pi)(pjPj) + (p^q) + ^S 2) {p\pj){p ] jPi) + {p-^q) 

+X i i ) e IJ e KL {p\pk)(p ] jPl) + (p -> q) + A 2 (P/P/)(9j?j) 
+^{p\tq*j){p T itqj) + *4,£ij£kl [(PjW)(Pj c 9l) + c.c] , (2.17) 



where c.c. denotes complex conjugate. Thus it is only Ai that proliferates — into three 
parameters — once flavor symmetric combinations are enumerated. In (12 ,4p one has 
the flavor specifications 



v\ (pUj) (qjPi) + v< i ] (pJ e Pj) (pWj) + {p^q) 

+ v {2) eue KL (pJep J )(p i K ep* L ) + (p -> q) 

+ V3eueKL[{pJepj){qKeq L ) + (p\ePj)(MK e 9l)], 



(2.18) 
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with a slight proliferation v% 



v, 



W „( 2 ) 



that is compensated by the fact that 1/4 = 



since that type of term always has a lone p or a lone q. In fact, this forbids the 
term for all cases Nf > 1. The 3/2 Yukawa term in ( 12.10[) is forbidden for the same 
reason as the z/4 potential term. The flavor symmetric Yukawa term is just (y± — > y): 

A/ = ?/[P/ A X P / + X/pAp/ + gjAxg/ + XgjAgi]- (2.19) 
Altogether one has in the Nf = 2 case the following twelve tunings to perform: 

• one scalar mass m 2 , 

. ten quartic couplings A?\ A?\ A< 3) , A 2 , A 3 , A 4 , ^, ,/«, ^, and 

• one Yukawa coupling y. 

As we will describe in more detail below, two approaches can be taken towards 
tuning y. In the first case, as was discussed above and advocated in [78], one adjusts 
the field strength of the scalars p, q in order to accomplish the same thing as fine- 
tuning y. In the second case, and this is a new approach that we propose for the first 
time here, one takes y to implicitly define the gauge coupling of the lattice theory 
and fine-tunes the Wilson gauge action coefficient (3 = 4/ ' g 1 until supersymmetry is 
achieved. The disadvantage of this second method is that one does not know a priori 
what the bare gauge coupling of the theory is really is! That is, it is determined in 
the process of offline reweighting. Yet since in a typical application all one really 
wants is say three values of j3 with sufficiently fine lattice spacing a, in order to make 
a continuum extrapolation, it reasonable to think that selecting three values of y will 
accomplish the same goal. In particular, one ought to bootstrap from small lattices 
where the actual value of j3 for the supersymmetry theory can be determined (in the 
reweighting process) cheaply. 



2.1.4 N f > 2 



The scalar mass terms are given by (12.161) . The quartic terms in (12 ,4p allow for seven 
flavor symmetric terms 

aS 1} (p\pi)(PjPj) + (p^q) + AS 2) (p\p.j)(p j jPi) + (p-*q) 

+AS 3) [t a utK L (p\Pj)(pyL) + (P- q)} + AS 4) [t a jjt a KL (p\p L )(pyj) + (P^ q) 

+Hp\pi)(q j jqj) + ^\pWj)(pJeqj) + X [ 3%t a KL (p\ep* K )(p T jep L ). (2.20) 
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Here t a u are generators of the flavor group SU(Nf). In (12.41) one has the flavor 
specifications 

vi(p\qj){<ijPi) + "2 (pJ e Pj){p\tp*j) + {p^q) 

+ [Vkl(pWk)(pWl) + (P- q)] ■ (2.21) 

The Yukawa couplings are given by (12.191) . Altogether one has twelve parameters to 
fine-tune: 

• one scalar mass m 2 , 

• ten quartic couplings A[ , . . . , \[ , A 2 , \ f\ Ui, , , an d 

• one Yukawa coupling y. 

2.2 SU(3) gauge theory 

Here triality (the Z 3 center symmetry of SU(3) gauge theory) is rather restrictive 
when combined with the other symmetries. The mass term is just as in the SU(2) 
theory above, Eq. (12.21) . Cubic potential terms are ruled out by £7(1) v. The quartic 
potential terms must be built from SU(3) singlet combinations of one of the formslfl 

fpjfpt, fPjt<lh <tw k <li- (2-22) 

because of U(l)y and triality; here ijkl are color indices (=1,2,3). To form color 
invariants we begin with enumerating the irreducible representations that occur from 
pairing: 

3x3 = 3 + 6, 3x3 = 1 + 8, (2.23) 

which in terms of fields takes the forms given in the Table [2j Some words of clarifi- 
cation are in order. First, conventional shorthands such as 

p ] p = p l Pil ph a p = p l (t a y iPj (2.24) 

and 

(pqT = 2 e * J *Pi9k> (pq)ij = 2 (PiQj + PjQi) (2-25) 



Here we follow the convention of raised indices on the 3 irreducible representation (irrep), and 
hence find it convenient to write p rather than p* 
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quadratic 


irrep 


Qv 


z 3 


p^p, q^q 


1 








q^p 


1 


2 





(pqY 


3 





2 


(ppY 


3 


2 


2 


{qqf 


3 


-2 


2 


(pq)ij 


6 





2 


(pph 


6 


2 


2 


(qq)ij 


6 


-2 


2 


pH a p, qH a q 


8 








qH a p 


8 


2 





pH a q 


8 


-2 






Table 2: SU(3) representations from pairs, flavor indices suppressed. 



have been employed. Here, {t a )l = (1/2) A?- with A a the Gell-Mann matrices. Second, 
flavor indices have been suppressed, but are necessary to render the compact notation 
sensible. For example, expressions such as e^ k pjPk require Nf > 1 in order to have 
nonvanishing result, ^ k Pj,jPk,K-, J,K = 1, . . . , Nf. Expressions such as piPj + pjPi 
would have a flavor specification PijPjj + Pj.iPi,j- For brevity, in Table [2] we have 
left out representations that can be obtained by complex conjugation, such as the 3 
representation e i jkp : 'p k - 

To obtain singlets we take the combinations which may be schematically denoted 
1-1, 3*3 i, 6 y '6y, 8 a 8 a . However two constraints relate these: 

{pq)i{pq) 1 = - [{p^p){q ] q) - {q j p){p^q)] , 

(pq) ij (pqh = \ Wp)(q ] q) + (q ] p)(p ] q)] ■ (2.26) 

Thus the 3*3j and 6 y- 6ij singlets can be eliminated in favor of the 1 • 1 forms, and one 
finds that the most general quartic Lagrangian is: 

A = \[{p ] p) 2 + (q ] q) 2 ] + Mp ] p)(q ] q) + Mp ] q)(q ] p) 

+H(ph a p) 2 + (qh a q) 2 ] + X 5 (ph a p)(qh a q) + X 6 (ph a q)(qH a p), (2.27) 
where as usual flavor specifications remain to be given (below), depending on the 
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value of Nt. Finally, the Yukawa couplings are 



C y = y[p T \Xp + Xp^P + Q T ^Xq + X T aM], 



(2.28) 



where A = \ a t a and A = X a t a , as in Eq. (12. lip . Note that the second type of term 
appearing in (12.101) is not allowed, due to SU(3) triality. 

2.2.1 N f = 1 



Here there is nothing to specify; the expressions just given suffice, with ( 12. 2p for the 
most general mass term, Eq. (12.281) for the Yukawa terms and (12.271) for the quartic 
terms. Altogether we have eight fine-tunings: 

• one scalar mass m 2 , 

• six quartic couplings Ai 5 2,...,6, and 

• one Yukawa coupling y. 

As before, the tuning of the Yukawa can be effectively accomplished either through 
the bare scalar field strength Z$ or through tuning the bare coupling (3 = 6/g 2 . 



2.2.2 N 



f 



A few operators proliferate because of different ways of realizing the SU (2) p x SU (2) q 
flavor symmetry. The quartic terms are: 



4 (p\pi)(PjPj) + (p -»• q) + a£ 2) (p\pj)(pjPi) + (p -> q) 
+\f ) t IJ t KL {p\pk){p ] jPl) + (p -* q) + M{p\pi){q ] jqj) 

+Hp\qj){q ] jPi) + $\(jp\rpiMfpj) + (p^ q)\ 

+\f[{p\t a pj){p\t a p I ) + (p^ q)} + \feijeKL\b\rpMrpL) + (p - ?)] 
+Hp\t a PiMt a qj) + K(p\t a qj)(q\t a pj). (2.29) 



The mass terms are given by (12.161) and the Yukawa terms by (12.191) . Altogether we 
have twelve fine-tunings: 



one scalar mass m 2 , 

(1 2 3) 

ten quartic couplings A2,3,5,6, A^ \ ' , and 
one Yukawa coupling y. 
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2.2.3 N f > 2 

We have in this case twelve operators in the quartic Lagrangian: 
£ 4 = \f> {p\pi){j)\pj) + (p -> q) {p\p.j){p ] jPi) + (p -> g) 

+Ai 3) [M £ (pJpj)(p!cPl) + (P - ff)] + Al 4) [M l (p}pz)(p!cPj) + (p - q) 

+Hv\pi){q\qj) + Hp\qj)(qjPi) + Ai 1) [(p}t»(p t J t> J ) + (p -> g)] 

+Af[(plt>j)(p t J t>/) + (p ^ g)] + AS 3) [« L (p^ a Pj)(pJ r t a p i ) + (p - g)] 
+Al 4) [t? J t^(p}t a PL)(p^ a Pj) + (p - g)] + Hp\t a P i)(qWqj) 
+Hp\t a qj)(q\t a Pi) (2-30) 



The mass terms are given by (12.161) and the Yukawa terms by f 1 2 . 1 9 j) . We now have 
fourteen fine-tunings: 



one scalar mass m 2 , 

twelve quartic couplings A^ , . . . , A] 4 , A2, A3, A 4 , 



j A4 , A5, A6, and 



• one Yukawa coupling y. 

Clearly the task of tuning these parameters such that the long distance effective theory 
is the much simpler Lagrangian flA.lj) poses an enormous challenge. A first task is to 
design a strategy for confirming from lattice data that the effective potential reduces 
to Eq. ( 12. 9ft . We leave this as a topic for future research. 



2.3 SU(4) gauge theory 

Here an analysis similar to what has just been performed for SU(3) leads to the 
conclusion that the only new quartic operator, not contained in (12.271) . is 



e l3kl piPjq k qi, + c.c. 



(2.31) 



In particular, using 4 x 4 = 6 + 10 to form quadratics in the 6 and 10 representations, 
one finds identities similar to (12.261) for 6 l - 7 6j 7 - and 10^10 3 that reduce these to 1 • 1 
forms. Thus the general quartic Lagrangian is: 

£4 = Ax[(p f p) 2 + (g f g) 2 ] + A 2 (p t p)(g t g) + A 3 (p t g)(g t p) 
+A 4 [(p t f I p) 2 + {qH a qf} + A 5 (p t t a p)(g t t a g) + \ & {pH a q){qH a p) 
+\ 7 [e i ' M p i p J q k q e + c.c.]. (2.32) 
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The mass and Yukawa terms are the same as for SU(3). 

The flavor specifications also follow SU(3). For Nf = 1,2, the counting of pa- 
rameters to be tuned is just increased by one relative to SU(3), due to the additional 
quartic coupling A7 in (12.321) above. For Nf = 2 it has the form 

h[tij£KL£ vM PnPjjqkKqeL + c.c.]. (2.33) 

However, that coupling is forbidden for Nf > 2 since it will not be SU(Nf) p x SU {Nf) q 
invariant. Thus for Nf > 3 the Lagrangian and parameter counting is identical to 
SU(3). 

2.4 SU(N > 4) gauge theory 

Here the form of the Lagrangian is just as in SU(3). The flavor specifications, de- 
pending on Nf are likewise identical. 

2.5 Summary 

In summary, SQCD contains O(10) fine-tunings in each case. For most of the the- 
ories, all of these tunings are bosonic and can be done offline. The exception was 
SU(2) with Nf = 1, where two Yukawa parameters must be adjusted. Setting aside 
that case, tuning between eight and fourteen couplings on bosonic operators will pose 
a significant challenge, even with the multicanonical reweighting techniques that we 
discuss below. A careful bootstrapping method, from small to large lattices, will be 
necessary in order to properly locate the critical parameter values in such a large 
parameter space. As mentioned above, it is best to begin with weak couplings on a 
small lattice, where lattice perturbation theory should be a useful guide. High statis- 
tics studies will be required in order to constrain such a large number of parameters. 
Tuning against the supersymmetry Ward identities, as will be described for the M = 4 
SYM case in the next section, also requires adjustment of mixing coefficients for bare 
operators appearing in the supercurrent. The problem appears daunting, and could 
only work if an automated, recursive simulate/search strategy is employed. 
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3 Lattice N=4 SYM 

In this section we describe a lattice formulation of four- dimensional J\f = 4 SYM0 
that may be within reach of practical simulations [78], when combined with the 
multicanonical methods described in §U As for SQCD, we use GW fermions to avoid 
gluino masses. Just as important, the GW fermions provide for an implementation of 
the global SU(A)n symmetry, which is chiral in how it couples fermions and scalars. 
The continuum chiral SU(4)r symmetry is replaced by a lattice generalization. As 
will be seen, this symmetry limits the number of counterterms that must be fine- 
tuned in important ways. As was the case in SQCD, only bosonic operators require 
fine-tuning; all tunings can be done "offline" by a Ferrenberg-Swendsen [81-83] type 
reweighting, exploiting multicanonical simulations to greatly broaden the parameter 
space that can be scanned offline. This aspect of the theory will be described in detail 
in §1 

3.1 Lattice Action 

The continuum field content corresponds to SU (N c ) Yang-Mills coupled to scalars and 
fermions in an SU(4) symmetric way. Typically one writes the global symmetry as 
SU(4)r, where the R denotes a symmetry that does not commute with the generators 
of supersymmetry. There are four massless Majorana fermions. The left-handed 
components transform in the fundamental 4 representation of SU(4)r and 6 real 
scalars in the 6 representation (antisymmetric tensor). If it were not for the Yukawa 
couplings, we could formulate the theory instead in terms of two Dirac fermions, 
which would simplify matters with respect to the GW formulation. However, the 
chiral ipRipipi ~ 4# • 6 • 4^ Yukawa couplings require that we decompose the fermions 
into four the left- and right-handed Majorana fermion components, which are related 
to each other by charge conjugation. (Note that ipR ~ 4 so that ipR ~ 4.) 

The six real scalars will be expressed with a single index (fi m , m—1. . .6, or composed 
into SU(A)ji Weyl matrices: 4>ij=4> m (Jm,ij and =4> m a i ^ l , where ex's are just SU{A) R 
Clebsch- Gordon coefficients involved in forming the singlet associated with 4# • 6 ■ 4^, 
or equivalently 43 6®4. They are most easily determined by dimensional reduction 
from ten dimensions, or by recognizing them as the six-dimensional Weyl matrices 
that are the building blocks of the six-dimensional Dirac gamma matrices. 

6 A review of the continuum theory, its superconformal representations and the AdS/CFT corre- 
spondence is given in [92]. 
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The Euclidean continuum action is 

S = 1 Tr | -G^G^ + {D^ m f + fapfr 

+V2$i (<f> ij P L -(<j> ij )*P R ) ^ + [(f> m ,<f>n)[<f>m,<f>n)j. (3.1) 

The SU(4:)r preserving scalar lattice action must allow for generic coefficients and 
non-supersymmetric terms, so that the supersymmetry-restoring counterterms can be 
tuned. The quartic interaction terms in the SU(2) and SU(3) case are: 

Ai Tr m 0„0 m n + A 2 Tr m m n n . (3.2) 

Comparing to (13.11) . we see that classically supersymmetry corresponds to 

Ax = 1/g 2 , A 2 = -1/g 2 . (3.3) 

In the case of SU(N C > 3), a total of four quartic terms should be included, both the 
operators (13.21) as well as 

A 3 Tr (j) m (j) n Tr m n + A 4 Tr m m Tr 0„0„. (3.4) 

For SU(2) and SU(3) these can be eliminated in favor of the single trace operators 
(13. 2\\ using algebraic identities. A scalar mass term must also be included: 

im 2 Tr 4> m (j) m . (3.5) 

Regarding the kinetic term (D M m ) 2 , one could use a naive gauge covariant nearest 
neighbor approximation. On the other hand, it has been seen in many previous studies 
that taking to be related to the fermion operator is advantageous to reducing 
supersymmetry-violating artifacts, presumably due to degeneracies of modes in the 
UV where weak coupling applies [54,93-96]. Obviously such an implementation would 
be more demanding numerically, since one uses the GW operator in the scalar sector. 
On the other hand, the advantages that might come in reducing lattice artifacts may 
well make it worth the effort. 

The precise type of Ginsparg- Wilson [55] fermion to be used, be they domain 
wall [97] or overlap [98], is not important for the considerations here. However, it 
has been argued that naive lattice Yukawa terms lead to inconsistencies in either the 
chiral or Majorana projections (depending on how the Yukawas are transcribed to the 
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lattice) [99]. Following Liischer [100], and Kikukawa and Suzuki [101], we introduce 
auxiliary fermions The fermionic lattice action is 



S F = 



^Tr {faDiPi-ViVi 



X 



+yV2$+*)i {(t> ij P L - {(t> ij )*P R ) ftH*);}- 



(3.6) 



where D is the GW operator. This action possesses an exact SU(A) R symmetry, with 
the scalars transforming as in the continuum and the fermions transforming according 
to 



operators, T is the generator of SU(4) R in the fundamental (4) and we have sup- 
pressed the SU(A) R indices. Hence (tp+ty) and (if+^f) transform like the continuum 
■0, fields. 

This auxiliary fermion method preserves the i?-symmetry exactly and keeps the 
Yukawa terms ultralocal. It is also consistent with the Majorana decomposition, so 
the fermionic determinant is an exact square; taking its square root to implement the 
Majorana nature of the fermions retains locality. The cost is the introduction of an 
extra fermionic excitation which is however nondynamical with 0(a~ v ) mass, so 
it decouples from the theory in the continuum limit. 

In the case of = 0, it is known that the overlap operator D is non-negative; 
in particular, det-D > 0. If the domain wall fermion approximation is used, then 
det-D > for this case. One can ask what happens to this positivity feature for 
7^ 0. It is easy to see that the fermion measure is real. In the field space (-0, the 
fermion matrix has the 2x2 block form: 



Sip/ie 

Sif/ie 
S^f/ie 



(TP L -T*P R )if), 

(T+T*) l5 DiP + (TP L -T*P R )^, 
(T*P L -TP R ) + ^(T+T*) 75 , 
(TP L -T*P R ) . 



(3.7) 



Here P L/R = |(1±7 5 ) 



^(1±75(1— 2D)) are the lattice modified chiral projection 




(3.8) 



Since 75.0*75 = D and similarly for My, we have 



(detAI)* = det.M f = det7 5 .M 1 7 5 = det M. 



(3.9) 
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The sign of the determinant may fluctuate. In fact, to agree with some results from 
the continuum (or really the zero- dimensional reduction — matrix models) we know 
that it must [102]. Thus a sign problem in the lattice theory may reflect continuum 
dynamics. One of the interesting questions in a lattice study is the extent to which 
this correlates with motion through the nontrivial moduli space in M = 4 SYM. In 
a simulation, the sign fluctuations will have to be accounted for by monitoring the 
low-lying eigenvalues of the fermion matrix, which can be computed efficiently 

3.2 Tuning to the supersymmetric theory 

Our goal is to nonperturbatively tune the lattice action such that the IR description 
is a good approximation to M = 4 SYM, with errors that are 0(a) and the lattice 
spacing a much smaller than the scales of interest. Due to operator mixing there is 
a nontrivial matching between the lattice and effective IR theories. All relevant and 
marginal terms consistent with lattice symmetries will appear in the infrared, except 
at special points in bare parameter space. We can arrive at the desired special point, 
M = 4 SYM, by introducing the supersymmetry-violating operators into the bare 
action and fine-tuning counterterms. These counterterms fall into three categories: a 
scalar mass term, a Yukawa term, and two or four scalar quartic terms, depending 
on the number of colors for the gauge group, restricted here to SU(N C ). As has 
already been mentioned, if N c < 3 then only two quartic terms need to be included, 
Eq. (13.21) . while for N c > 3 two more quartic terms must be introduced, Eq. (13. 4p . 
As in the SQCD discussion above, rescaling the Yukawa term can be accomplished 
through a rescaling of the scalar kinetic term. Therefore in the lattice theory the 
scalar kinetic term should be taken to have a general coefficient that is to be 
tuned nonperturbatively. We will describe the multicanonical reweighting method 
of fine-tuning the M = 4 SYM lattice theory in more detail in §H below. This will 
include a discussion of mixing coefficients that must be measured in the supercurrent 
in order to use supersymmetry Ward identities in the fine-tuning process. 

We now comment on the significance of preserving the chiral symmetries of the 
theory, albeit in the lattice-modified form (13. 7p . For this purpose suppose we were to 
formulate the fermionic part of the theory instead as 

s F = Tr + 4^ (4> ij PL - {</> iS yp R ) ^} (3.io) 

X * 

with D w the Wilson-Dirac operator. Then due to the explicit violation of chiral 
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symmetry, a mass correction 

rmpiipi = m (ipu^Ri + 4>m4>Li) (3.11) 

would be generated. Note that the Majorana condition ip^ = ipi implies that ipu and 
ipRi are in conjugate representations: 

P R r = C$1 = ipRi P L if> c = C^l = ipL, (3.12) 

Comparing to (13.111) . we see that the mass term consists of 4 • 4 and 4 • 4 couplings, 
so that only the real subgroup SO (A) of the original SU{4)r symmetry is preserved 
in this Wilson-Dirac formulation. Futhermore, this mass term converts ipu into ipm 
(thus mixing the 4 and 4 representations of the SU(4)r), so we will radiatively gen- 
erate new Yukawa couplings 

y'yftfa (<j> ij P R - (<t> ij yP L ) V;, (3.13) 
where chiralities have been swapped relative to the supersymmetric Yukawa term in 

On one hand, the SO (4) that is preserved does limit the number of parameters 
that must be fine-tuned. We have just two more, m and y', than in the GW case. 
On the other hand, these are additional fermionic counterterms, and we cannot use 
the trick of rescaling the scalar kinetic term, since that freedom has already been 
exploited for tuning the parameter y. So, we face the problem that two new terms of 
a very problematic type are present because we have not preserved the chiral SU (4)^, 
but only a real subgroup. This will also pose a problem for tuning of Ward identities, 
because more operators can mix with the supercurrent, due to the reduced symmetry 
of the lattice theory. That translates into additional mixing coefficients that must be 
measured nonperturbatively. 

3.3 Discussion 

We have used GW fermions with chiral SU(4) R invariant Yukawa couplings following 
the method of [100, 101] to reduce the counterterms in an important way. Because 
the counterterm tuning bosonic, it can be done offline. In §4] we will explain how to 
alleviate the ensemble overlap problemby taking a multicanonical approach, flattening 
the distribution with respect to the parameters that are to be scanned. 

We invite the reader to contemplate the following circumstance, which is interest- 
ing since it is radically different from what occurs in lattice QCD: since M = 4 SYM 
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is conformal, the continuum limit is not a weak coupling limit g — > 0. What is known 
about such lattice field theories? Perhaps the best starting point is the class of two- 
dimensional models with an IR fixed point that have been extensively on the lattice. 
But M = 4 SYM should have not merely an IR fixed point (as is supposed to occur 
in some 4d gauge theories that have been studied recently on the lattice [103-105]), 
but a whole critical line of fixed points. The well-known two-dimensional analogy is 
the XY model (planar spin model, 0(2) model, etc.). In that case, each k > k c leads 
to a conformal field theory (CFT) in the IR, but it is in fact a continuous family of 
CFTs, with scaling exponents (anomalous dimensions) that depend on kJ^| The same 
should be true for M = 4 SYM: each value of the continuum gauge coupling g CO nt. 
corresponds to a CFT, but the anomalous dimensions of non-chiral, or non-BPS, 
quantities will depend on the value of g CO nt.- The lattice formulation should also have 
that feature, though the relation between the lattice coupling gi att and the continuum 
one g con t. must be established through detailed calculations. 

Our proposal should certainly work at very weak coupling gi aU _ ^ g CO nt. "C 1, where 
one knows that the IR description will be in terms of the same degrees of freedom 
as one puts on the lattice. Here it is important that, due to the fine-tuning, one 
starts with gi att _ at the lattice scale and the theory flows into the IR fixed point such 
that the coupling remains weak, terminating with the value g CO nt. < 1. Of course 
in such a situation perturbation theory is reliable and there is no need to use lattice 
discretization. Nevertheless, this is a useful reference point for the lattice parameters 
that must be fine-tuned. 

But is it guaranteed that one can find lattice parameters which correspond to the 
strongly coupled continuum theory? No definitive answer to this question will be 
offered here, though we certainly have our own opinion. We only see room for four 
possibilities: 

1. No lattice field theory can describe strong TV = 4 SYM, but there is another, 
more sophisticated nonperturbative formulation that can do the job. 

2. Strong Af = 4 SYM does not really exist — it is just a continuum field theorists' 
fantasy. 

3. A lattice field theory exists that will do the job of describing the strongly coupled 

7 Recall that large k corresponds to low temperature, which is precisely where order is anticipated. 
Of course in the XY model it is just algebraic order, due to the absence of spontaneous symmetry 
breaking in two dimensions. 
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theory, but it is not the one formulated here. 



4. Along a line of points in the parameter space of the lattice theory described 
here, strongly coupled M = 4 SYM emerges in the IR. 

Our opinion is that first two possibilities are virtually impossible. In particular, 
there is convincing evidence from the AdS / CFT correspondence that strong M = 4 
SYM corresponds to a weakly coupled supergravity theory, and that it is perfectly 
consistent [45-47]. In the gauge theory, the anomalous dimensions of BPS operators 
are determined by the superconformal algebra and are protected from renormaliza- 
tion. Thus they can be computed at arbitrarily weak coupling and then continued 
into the strong coupling regime. Thus certain operators written in terms of elemen- 
tary fields have a well understood behavior in the IR. Given that this is true, it is 
hard to imagine why the elementary fields would not provide "the correct degrees of 
freedom" in terms of which to define the theory at strong coupling. 

Furthermore, the renormalization group perspective indicates that some lattice 
theory should exist whereby the lattice artifacts can be cancelled by irrelevant oper- 
ators, leading to a perfect lattice actional Such a perfect lattice action theory would 
only differ from the one we are proposing by (naively) irrelevant terms. The fine- 
tuning approach that we are advocating could only fail (Option 3 in the list above) 
if the naively irrelevant operators that appear in the perfect lattice action theory 
turn out to be relevant or marginal in the IR as one moves to stronger lattice gauge 
couplings. 

However, we believe that Item 4 is the true state of affairs, though at this point it 
is a matter of speculation. The lattice will introduce artifacts that break conformal 
symmetry by irrelevant operators. This will be characterized by the lattice spacing 
a. Thus starting at the lattice scale, there is a renormalization group flow into the 
IR fixed point where the theory becomes conformal. Equivalently, one must look at 
distances i > a in order to see the conformal field theory behavior. 

For instance, correlation functions over a distance x will depend on a and the 
linear size of the lattice L, where the latter serves as an IR regulator. They will have 
the general form 



8 Of course this assumes that continuum power-counting can be applied. However, it is reasonable 
to suppose that with sufficient effort one's intuition in this regard can be proven, as has been done 
recently for staggered fermions [106]. 




(3.14) 
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where h is the scaling dimension of the operator and the 0(a/x,x/L) correction 
carries the scaling violation due to irrelevant lattice artifacts, as well as the finite size 
corrections. One cannot simply take L — > oo to remove these, since the theory needs 
an IR regulator if it is to make sense, under the assumption that it flows into an IR 
fixed point. 

Thus, suppose the properly tuned lattice action flows into an IR fixed point. There 
is some scale ^cft beyond which further flow is negligible. For x > £qft, the behavior 
is that of a CFT. In order to capture this regime, it is important that £qft 

< L. We 

know that we have a fine lattice (small lattice spacing a) if 

£ C F T /a > 1. (3.15) 

To identify the distance scale at which conformality sets in, one could monitor the 
running gauge coupling using Schrodinger functional and step-scaling methods [107, 
108], as has been done in [103-105, 109]. 

To have a strong CFT, one should look for 0(1) anomalous dimensions through 
measurements of critical exponents. Admittedly this would be an enormous challenge 
to accomplish through the standard method of finite-size scaling studies. This is 
because the theory is four-dimensional with GW fermions and fine-tuning that grows 
exponentially more difficult as the volume is increased. On the other hand, one can 
fantasize that it may be possible to extract information from small volume studies, 
analogous to the e regime methods that are used in QCD. If one were to accomplish 
this super-human task, find 0(1) anomalous dimensions, and verify that all of the 
£77(2, 2|4) x SU(A) R Ward identities are satisfied, then it is hard to see how one 
has anything other than strong M = 4 SYM, based on universality arguments. The 
capstone of such a study would be to match the anomalous dimensions to those 
predicted by the AdS/CFT correspondence. 

In general, lattice theories have a phase structure that is richer than the continuum 
theory that they are intended to define. A well-known example is SU(2) gauge theory 
with a mixed fundamental/adjoint Wilson action [110,111]. Phase boundaries exist 
and some regions of the lattice parameter space have only 0(a) correlations. These are 
"lattice phases" separated from the phase with a continuum limit by a bulk transition. 
This will happen in the M = 4 SYM SU(2) lattice theory, since several adjoint 
fermions are present; under coarse- graining they will generate adjoint Wilson action 
terms: 

S eff 3-(3 A ^^U^\x). (3.16) 

X 



24 



Indeed, this has recently been observed to give rise to a bulk transition in SU(2) 
lattice gauge theory with four adjoint Majorana fermions [112]. Similar findings have 
been reported for SU(3) with sextet fermions [113,114]. The situation in M = 4 SYM 
lattice theory with gauge group SU(2) will be similar to what was found in [112]: a 
continuum phase exists only if the bare lattice coupling (3 = 4/g 2 is sufficiently weak 
(e.g., g 2 < 2 in [112]). 

The point here is that even though the target theory is one in which the coupling 
constant does not run, in the lattice theory irrelevant operators cause important 
renormalizations at short distances so that the true strength of the IR coupling must 
be determined by a detailed study of the long distance physics. 

4 Fine-tuning with multicanonical reweighting 

Multicanonical methods [84-86] combined with "Ferrenberg-Swendsen reweighting" 
[81-83] [refered to here as multicanonical reweighting (MCRW)] have proven to be 
a powerful tool for maximizing the usefulness of Monte Carlo simulations over a 
range of parameter space much wider than was actually simulated. For instance, 
MCRW was applied in a study comparing SU{2) and 5*0(3) = SU(2)/Z2 lattice 
gauge theories [88,89]. It was found to dramatically flatten the distributions with 
respect to three parameters, twists on gauge fields at the spatial boundaries. Another 
successful application of MCRW consists of lattice results for the electroweak phase 
transition [87, 115]. 

We will begin by describing MCRW generally, followed by a presentation of how 
it would be applied to the Af = 4 lattice SYM that was described in the previous 
section. 

4.1 Preliminaries 

Suppose we perform a Monte Carlo simulation at one value m of the scalar mass to, 
so that the configurations sample the distribution determined by the action 

S(m = mo) = S(m = 0) + ^ J d 4 x m 2 <f) 2 (x). (4.1) 

Following the "Ferrenberg-Swendsen reweighting" method [81-83] one can use the 
following reweighting identity to compute the expectation value (0) m of an operator 
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O for the distribution with a mass m: 

(n) _ (Oe- As ^) mo Ec ef( n)Ocexp ^-rjjfe 

In the first equality (■ ■ •) is the expectation value with respect to the canonical 
distribution corresponding to (14.10 and 

A5(m) = i(m 2 - m 2 ,) J d 4 x 2 (4.3) 

is the shift in the action when the mass is changed. In the second, J d 4 x (p 2 c and 
Oc are the mass term and operator evaluated on configuration C and Ec*eF(n) ^ s 
the sum over the distribution F(n) of n configurations generated in the Monte Carlo 
simulation. These of course provide a finite ensemble that approximates the canonical 
distribution corresponding to (14.11) . The advantage of this approach is that one need 
only perform a single simulation at mass m , storing the values of J d A x 4> c and Oc 
for each C, and then (O) can be computed for a swath of the parameter space m 
without having to perform any new simulations. Typically the time for this "offline" 
calculation is negligible compared to that of the simulation. 

Unfortunately, the regime of utility for this technique is limited by the overlap 
problem, in a way that often worsens exponentially in the spacetime volume. For 
instance, suppose the theory (14. ip has a quartic interaction and a critical mass-squared 
m 2 such that for m 2 < m 2 c there is spontaneous symmetry breaking. If we simulate 
with rriQ > m 2 then the field is exponentially weighted toward J d 4 x <ft 2 ~ 0. Now 
suppose we attempt to reweight to m 2 < m 2 . In that case — (m 2 — wip) > so that 
the exponential weight factor in (14.21) is minimal at j d 4 x <fi 2 ~ 0. The ensemble that 
is generated in the Monte Carlo simulation will have exponentially few configuration 
in the regime where J d 4 x <p 2 is far from zero and e -A5 ( m ) is large. Because we will 
have very few representatives of configurations with the largest weig ht e" AS(m) , and 
most members of the ensemble have very small weight, fluctuations will be large and 
huge samples are required in order to have acceptable errors. The mismatch of the 
distributions gets worse as the number of lattice sites increases, because the exponent 
is extensive (i.e., scales like the spacetime volume L 4 ). 

As a concrete example, we reproduce in Fig. [1] a figure from [89] . It shows that in 
the range of a three-dimensional parameter space the ordinary canonical Monte Carlo 
distribution varies by 15 orders of magnitude. This is for an 8 3 x 4 lattice, which is 
still relatively small. 
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Figure 1: Fig. 4 of [89], with permission. This shows the unweighted distribution (blue 
burst), multicanonical reweighting function (red plus) and the reweighted distribution 



sreen x . 



In a number of contexts the technique of multicanonical reweighting [84-86] has 
been found to ameliorate the overlap problem. One replaces S with 

S MCR w = S + W[0 1 ,0 2 ,...}, (4.4) 

where W[Oi, O2, ■ ■ ■] is a carefully engineered function of some small set of observ- 
ables. For instance in the M = 4 SYM case W will be a function of J d 4 x <p 2 , the 
distinct quartic terms J d A x 4 and the kinetic term j d 4 x (D<p) 2 . The (reweighted) 
expectation value of an observable in the distribution corresponding to Smcrw is: 

(O) = ^ C&F{n) ° (4 5) 

Since the e w factor in (14. 5p just cancels the e~ w Boltzmann factor coming from 
f)4.4p . one might wonder why it is introduced in the first place. The point is that 
the additional Boltzmann factor e~ w in effect produces a weighted average over a 
continuum of canonical ensembles (hence the appelation "multicanonical" ) such that 
there is a good overlap with the distribution that one is reweighting to. The challenge 
is to design a W such that sampling is flattened over the range of observables one is 
interested in. 
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We return to Fig. [TJ taken from [89]. It shows that the multicanonical Monte 
Carlo sampling distribution is flat in the range of three-dimensional parameter space 
between the peaks, where the ordinary canonical Monte Carlo distribution varies 
by twelve orders of magnitude. The reweighting function W was represented by a 
numerical table, composed of the inverse density of states with respect to the tuned 
parameters. This is for an 8 3 x 4 lattice, which is still relatively small, and it indicates 
that O(10 12 ) more samples would be required in the canonical Monte Carlo approach 
in order to scan a comparable range of parameter space by ordinary Ferrenberg- 
Swendsen reweighting techniques. Working on lattices of size, say, 16 4 , would make 
the overlap problems of the canonical distribution many orders of magnitude worse. 
But lattices of this size and larger are needed in order to extract continuum behavior 
from the lattice. On the other hand, it is not known how difficult the overlap problem 
is in the two types of supersymmetric models considered above, SQCD and M = 4 
SYM. 

As another example, in studying first order phase transitions (e.g., [87]), one 
chooses 0\ to be the order parameter of the transition; in a model with a scalar field, 
typically 0\ = J d 4 xcf) 2 . One tunes VK[C?i] to cancel the nonperturbative effective 
potential for this operator, so that the Monte Carlo simulation samples evenly in 0\. 
This enhances statistics for configurations intermediate between the phases. In the 
mass scan example of Eq. (14.21) , one has 

EceFM Oc exp (W[J d*x <%\ -\{m\- m\) J d*x 
{ } Ecmn) exp (W[J d±x $1] -\{m\- m\) J d*x <f> 2 c ) 1 " J 

In this way, wherever the exponential in (14.61) happens to be at its maximum, a large 
number of configurations will be generated, due to the flat distribution with respect 
to J d 4 x (j) 2 . 

Two approaches exist for engineering a good function W. 

(1) One can employ a bootstrap method that iterates between Monte Carlo simu- 
lation and adjustments to W . For instance a numerical tabulation of density of 
states p may be obtained from a canonical simulation, as was done in [88,89]. 
Schematically, one obtains a histogram estimate of p{0\) for an operator value 
range (D\ range 0i )Tn in < 0\ < Ox jmax . This provides an initial version of W, 
through W = 1/ p{0\). If necessary, the process can be repeated to refine the 
table. 
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(2) Iterative or stochastic searches may be used to optimize W with respect to 
a predetermined parameterization in a small volume. Performing this at two 
different small volumes then provides an extrapolation estimate for W in the 
next largest volume, which can then be refined through another search. 



4.2 Application of MCRW to M = 4 SYM 

Numerical studies of MCRW for M = 4 SYM have yet to be attempted, though the 
groundwork for this effort was laid in [78]. Here we will review those findings, as an 
illustration of the MCRW approach to lattice supersymmetry. 

For N c < 4, the reweighting function W will depend on the four bosonic contri- 
butions to the action 

Oi = d 4 x Tr D^mD^m, 2 = d 4 x Tr m m , 



C>3 = / d 4 X Tr <p m <Pn<Pm<Pm A = J d^X Tr <\) m <\>m<\>n<\>n- (4.7) 

For N c > 4 we must also include the double-trace operators 

05 = J d 4 x Tr m n Tr (j) m (j) n , C 6 = J d A x Tr m m Tr n n . (4.8) 
A finite sampling of the multicanonical ensemble described by 

Z MC = j [dA # d<f>) exp [-S(Zf\ ml Af) - W^O,})] (4.9) 

is performed by Monte Carlo simulation, where Z^\ m^, \f*' are the coefficients of 
the operators (14. 71) and (14.81) that appear in the action. In the simulation the RHMC 
algorithm [116-118] should be used as this is currently the best available approach 
for dynamical fermions. Of course one replaces the continuum gauge fields A^(x) 
with lattice link fields U^x), and the fermions are replaced with pseudofermions 
X l ( x ), with a corresponding reformulation of the lattice action, in the usual way. One 
then computes reweighted expectation values of quantities O for a different set of 
parameters Z^, m 2 , Aj via the relation 



(O) 



(Oex V [-AS + W({O t })]) MC 
mcrw {exp [-AS + W({O t })]} MC ' 



AS(Z , m\ \i\{Oi}) = S{Zt, m\ \\A, V, <p) - S(Z™ ,n%, AHA V>, 0) (4-10) 
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where MC indicates the multicanonical expectation value following from simulation 
with (14. 9J) . We have made it explicit that the change in action AS* depends only on 
the operators Oi that were given in (14. 7p and (|4.8p . 

The finite sample F(n) generated from the Monte Carlo simulation with distri- 
bution described by (14.91) consists of a set of configurations C\, . . . , C n . Thus for a 
given C the lattice fields U^x), X l (x), <fi m (x) take values U fl (C\x),x l {C\x),(f) m (C\x). 
From this point of view the MCRW evaluation of expectations values (I4.10p can be 
interpreted as computations with partition function 

^ MCRW — lim > exp[-AS c ({C?i}) +W c ({O i })} (4.11) 

ra— >oo ' ' 

CeF(n) 

where the subscripts C inside the exponential indicate that all fields are to be evalu- 
ated on configuration C . 

Next we define a density of states n(Oi, . . . , O4) = n({Oi}) in the multicanonical 
distribution, where we have specialized to the N c < 4 case for notational simplicity. 
Let $({Oi}) be any function of the operators. Then: 

J [dA # d<f>] e s[A^]-w[{o m] ^{o^]}) = 

J [dA # d4>] e -siA^]-w [{ o m] ( J YldO^SiO^] - Cf ) 

fldOf }n({OF})$({0®[<l>]}) (4.12) 
Thus we can write the reweighted multicanonical partition function (14.111) as: 

Zmcrw — 

w(Z (f) ,m 2 ,X i \{0, l })=ex P [-AS(Z (t) ,m 2 ,X i \{O i }) + W({O l })] (4.13) 

The engineering of iy({Oj}) has as its goal the generation of ensembles F{n) such 
that there is a reasonable number of configurations n({C?f }) with large weight for 
a broad patch in the parameter space Z^,m 2 , A, that we intend to scan in the fine- 
tuning process. For a choice of Z^, m 2 , A« within that patch, what we therefore want 
is n({Cf°' ) }) not too small wherever w(Z^,m 2 , Xi\{Oi}) has most of its support. Of 
course one must also decide where the patches of interest lie. This is best achieved 
through the bootstrap method described above (i.e., starting from small lattices and 
weak couplings, where the counterterms can be determined reliably using analytic 
methods). 
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4.2.1 Effective potential 

Tuning the scalar mass term. We begin by considering the case where we only 
shift r?r.Q — > m 2 relative to the reference point {Z^\ ttIq, Xf) of the multicanonical 
ensemble. One sees from f |4. lOj) that AS* = (m 2 — mfyOz, where O2 is the mass 
operator defined in f)4.7p . The gauge invariant effective potential in finite volume is 
defined as follows: 

e -nv cS (m^) = ^ e w({c?0)-(m 2 -mg)O 2 6 ^ A 2 _ &2\ \ ^ M ) 

where Q is the spacetime volume. Thus A 2 represents the mean value of the squared 
scalar field Tr (f> m (f) m . 
Now suppose we vary 

m 2 ^ m 2 + Am 2 . (4.15) 
In ( I4.14p we can use the 5-function as follows: 

g-Am'ftj f A 2 -Q\= e ~ AM2 5 (a 2 - ^ . (4.16) 

Since we can take e - Am2nA2 outside the expectation value, it is clear that the shift 
(14.151) changes V e g(m 2 \A 2 ) by adding a linear component: 

V cS (m 2 + Am 2 1 A 2 ) = V cS (m 2 \A 2 ) + Am 2 A 2 . (4.17) 

Therefore, measuring V e f[(m 2 \A 2 ) immediately determines how (Tr 4> m (p m ) and F M qrw 
— In Zmcrw vary as a function of m 2 . 

For example, suppose we measure V e s with m 2 = 0, and obtain a parameterization 

V eS (m 2 = 0\A 2 ) = rA 2 + uA 4 + ■ ■ ■ (4.18) 

where • • • represents terms higher order in A 2 . Then it is clear that the critical mass 
is m 2 = —r. 

As an alternative, one can also locate the critical m 2 by looking for the peak in 
the 4> m susceptibility 

X {4>) = J d A x (Tr m (x)0 m (O)> (4.19) 

On the lattice, gauge-fixing will be necessary since the fields are located at different 
sites x. As with the determination of the effective potential, this can be done offline 
using the reweighting techniques. The peak of the susceptibility x{4 > ) with respect 
to m 2 should agree with the point located by the effective potential analysis that we 
discussed in relation to (14.141) . 
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Tuning the quartic terms. Finding the quartic parameters that lead to additional 
second-order behavior should be possible. Indeed, it has been successfully achieved in 
the context of the electroweak phase transition [87]. In the M = 4 theory one could 
for instance look for peaks in the following susceptibility tensor: 



Xmn;pq J d X {O mn (x) O mn (0) ) conn j 

O mn = Tr (p m (pn - ~Smn ^ Tr (p k (p k , (4.20) 

k 

where "conn." denotes the connected correlation function. On the basis of SO(6) and 
cyclic trace symmetries, one has for the susceptibility: 

X,mn;pq C\(5 m p5 n q ~\~ $mq&rip) ^2^mn^pq- (4 .21) 

In the supersymmetric theory, the operator in (14.201) is the chiral primary operator 
of the supergravity multiplet, where the terminology arises from the AdS / CFT corre- 
spondence. It is 1/2 BPS so that its conformal dimension A = 2 is protected, and 
the susceptibility must reflect the fact that it transforms as a 20' of SO (6). Taking 
these properties into account, we have at the supersymmetric point the prediction 



^mn\pq 



{8 m p8nq + ^mq^np) hl(L/ d) , (4.22) 



where L is the linear size of the lattice, which serves as an IR cutoff. The term C2 in 
(I4.2ip vanishes in the supersymmetric limit since O mn creates an exact chiral primary 
state. Of course we work here with bare lattice operators and mixing will occur. 
However, the logarithmically divergent susceptibility should provide a clear signal of 
the additional second order behavior associated with tuning to the supersymmetric 
point. It is interesting to contrast this with the scalar susceptibility of the Konishi 
multiplet, Ok = J2k ^ r 0fc0fc, which according to perturbation theory [119] has Ak > 
2, and hence finite susceptibility with respect to L — > oo at fixed lattice spacing 
a. Because it is non-BPS, the AdS/CFT correspondence cannot be used to check 
the conformal dimension Ak at strong coupling. This is one feature that might be 
probed with the lattice, by examining the finite-size scaling of the corresponding 
susceptibility. 

Thus suppose that as we tune the quartic parameters Aj toward the supersym- 
metric point, flat directions of the scalar potential open up, revealing the moduli 
space of Af = 4 through divergent scalar susceptibilities. This presents numerical 
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difficulties and one would rather see the divergence as a limiting behavior. While it is 
true that lattice artifacts regulate this divergence, it is nevertheless desireable to have 
an independent knob that controls it. Furthermore, near the super symmetric point 
the quartic potential can turn over, leading to a runaway instability. To regulate the 
runaway directions and render susceptibilities finite, we propose to add a sextic term 



to the potential. Here, a is the lattice spacing, or in a continuum description, the 
inverse of the UV momentum cutoff. One might worry that radiative corrections 
in the lattice theory could cancel this term so that instabilities would not be cured. 
However, the instabilities are associated with large scalar field values, where the gauge 
symmetry is effectively broken. Hence the runaway directions correspond to sectors 
of the theory with weak gauge coupling, and so if the sextic coupling /i 6 appearing in 
(14.231) is sufficiently large the radiative corrections cannot overpower the stabilizing 
term (14T23D . 

In fact, the flat directions imply that unbroken Af = 4 SYM is not well behaved 
in the continuum limit at finite volume; the moduli are not fixed and the partition 
function diverges because of the integral over the infinite moduli space. Therefore it 
will always be necessary to break supersymmetry somehow. In addition to adding 
(I4.23p . we advocate introducing antiperiodic boundary conditions for the fermions 
in the temporal direction — finite temperature. This lifts the moduli degeneracy in 
a way that is removed by a zero temperature extrapolation, such as working on an 
L 3 x 2L lattice, with 2L/a sites in the temporal direction, and scaling L — > oo. The 
antiperiodic boundary conditions will also be beneficial to numerical stability of the 
dynamical fermion algorithms. 

4.2.2 Tuning with supersymmetric Ward identities 

If supersymmetry is exact then the supercurrent S^i is conserved. The index i cor- 
responds to i?-symmetry, with the supercurrent transforming as a 4 with respect to 
this group. It follows from this supercurrent conservation law that (d ll S fl ^(x)0(y)) 
vanishes at x^y for all local operators O. We can use this property to fine-tune to 
a supersymmetric point in parameter space, a technique pioneered in Af = 1 super- 
symmetry with Wilson fermions by the DESY-Miinster group [120]; here we discuss 
the extension to Af = 4 SYM. 




(4.23) 
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In the continuum, the supercurrent S^j is a linear combination of three dimension- 
7/2 operators. It is easy to write down corresponding lattice operators, though they 
not unique. In fact since we work with a cutoff theory, the lattice operators mix with 
continuum operators of higher dimension in the same symmetry channel. For this 
reason we express the operators O^i in a continuum language, since our purpose is 
to convey the method rather than the details of a specific implementation. 

An analogous mixing analysis occurs in Wilson fermion lattice Af—1 SYM. The 
DESY-Miinster group found that two dimension-7/2 operators, the supercurrent 
and another fermionic current T M , mix in the lattice-continuum matching. For this 
reason it is necessary to combine two corresponding lattice operators with unde- 
termined coefficients in order to find the lattice operator that becomes in the 
continuum limit. 

In our case we found in [78] that five dimension-7/2 operators must be taken 
into account. We denote them as Oj^'- 5 , and the renormalized supercurrent is in all 
generality of the form 



The terms on the right-hand side are bare operators. At tree level the supercurrent 
corresponds to Zi=Z 2 =Z 3 =l and Z^=Z 5 =0. The renormalization constants Z n are 
universal with respect to the index i due to S77(4).r symmetry. 

We wish to find the point in parameter space where the renormalized current of 
(I4.24p satisfies d^S^i = as an operator relation. One therefore measures correlation 
functions containing d^S^. In actuality, we demand supercurrent conservation up to 
0(a) corrections, since at finite lattice spacing there will always be supersymmetry- 
breaking due to lattice artifacts. What we seek is a trajectory in parameter space 
such that supersymmetry-breaking that vanishes in the continuum limit. 

Consider the SU(2) case, and suppose we have already tuned the bare mass m 2 
and the ratio of quartic couplings A 2 /Ai using the effective potential and susceptibility 
methods described in §4.2. II above. This leaves two more fine-tunings Z^ and Ai to be 
performed using the supersymmetric Ward identities. To tune these two parameters 
we need to examine correlation functions of six operators 0" i in the same symmetry 
channel as S^. The natural choice is the set C^'i ' 5 appearing in (14.241) . plus one 



-z* I** + m%Pi) }l^A 

+ {Z Alu F^5 i0 - Z 5 V2D^P L + fajPn^A + 0(a) 




(4.24) 
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dimension-9/2 operator & ^ One then measures the matrix of correlation functions 

M mn (t) = J rf 3 x(O m /(t,x)^(0,0)) (4.25) 

The t derivative of this is the correlation function between d^O^ and Oq i at vanishing 
spatial momentum. 

For the dimension-7/2 operators appearing in (14.251) . correlation functions will 
generically decay as 

(<(t,x)(^(0,O)) ~ (t 2 + x 2 )- 7 / 2 . (4.26) 

Integrating with respect to x, we therefore find that these elements of M mn (t) to 
decay as 

M mn {t) ~ r 4 . (4.27) 

At the supersymmetric point and for the right choices of the coefficients Z m appearing 
in (I4.24p . the corresponding combination of correlation functions Z m M mn will be 
suppressed by the lattice spacing a and hence decay as at~ 5 for the n associated with 
dimension-7/2 operators. In fact, GW fermions are automatically 0(a) improved, 
so if the operators were likewise improved we could even achieve a suppression a 2 t~ 6 
at the supersymmetric point, which would be easier to distinguish from the generic 
t~ A behavior. Given the cost of the GW simulations, and the fact that operator 
improvement is performed offline, it would be well worth the effort. For the dimension- 
9/2 operator we would require decays with at~ 6 , or a 2 t~ 7 if improvement is performed. 

These six conditions on the correlation matrix M mn (t) fix the four ratios Z2...5/Z1 
and the two parameters Z^, Ai that must be fine-tuned for the SU(2) and SU(3) 
cases. In the SU(iV c > 3) case, additional tunings with correlation functions of the 
supercurrent will be necessary. 

4.2.3 Other Ward identities 

In the superconformal phase of M = 4 SYM, (<p m ) = 0, the global symmetry of 
the theory is the supergroup Si7(2, 2|4) and the R-symmetry group SU(A)r. The 
lattice preserves the latter in a modified form, but deviations from the local (contin- 
uum) form could perhaps serve as a measure of the lattice spacing. The supergroup 
SU(2, 2|4) includes conformal supercharges 0.% other than the four supercharges Qi 
corresponding to the supercurrents S^i in (14. 24ft . The four additional supercurrents 
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@ M i could also be used as a probe for the M = 4 SYM theory, since they will have their 
own Ward identities. On the other hand, conservation of combined with scale in- 
variance and Poincare symmetry implies conservation of (5^, so measurements of the 
additional Ward identities are not independent, but rather serve as a means to check 
consistency with predictions of the continuum theory. Verification of this feature 
would be reassuring in the regime of strong IR gauge coupling. 

4.2.4 Summary 

We have seen that for SU(2) and SU(3), there are four fine-tunings in the action, 
Z^,m 2 ,Ai,A 2 . For SU(iV c > 3) colors there are six, Z^, m 2 , Ai, . . . , A 4 . In addition, 
one must fix the four relative renormalization constants Z 2 /Zi, . . . , Z 5 /Zi in the su- 
per current. 

It is conceivable that all but one of the scalar potential counterterms Aj can 
be fixed by matching the effective potential to the supersymmetric scalar potential 
Tr [<f) m , <f) n ] [cj) m , 4> n ], though the practicality of this is yet to be established. The overall 
strength of the quartic terms cannot be determined from the effective potential. The 
critical m? will be determined from the effective potential. 

Thus we see that in the more optimistic scenario, where the effective potential 
can be fully exploited, only two parameters of the lattice action need to be fine-tuned 
by the supersymmetric Ward identities: one fine-tuning of the bare kinetic coefficient 
Z^ for the scalar, one overall scalar potential coefficient Ai, and the four relative 
supercurrent coefficients Zi/Z 2) ■ ■ ■ , Z 4 /Z 5 . Hence a total of six Ward identities must 
be measured. 

If it proves too difficult to constrain all the ratios Aj+i/Aj of the quartic terms 
using effective potential methods, then additional supersymmetric Ward identities 
must be measured. An intelligent strategy would be to perform a combined mini- 
mization of all quantities that can be measured with reasonable accuracy, so that in 
fact the adjustment of parameters and mixing coefficients is overconstrained. Here, 
the additional Ward identities mentioned in §4.2.31 could also be employed. 

Aside from the challenges of developing an effective, automated and optimized 
tuning strategy, one must also face the fact that the lattice formulation employs GW 
fermions, which are numerically expensive. This is particullary true in a theory such 
as this one, with massless fermions and the corresponding critical slowing down. 

A first computation that needs to be done is to fix the multicanonical reweighting 
function. Small lattices (4 4 , 6 4 ) should suffice to get a rough idea of how to proceed in 
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further studies (8 4 , 10 4 , 12 4 ). Perturbative calculations may also help to narrow the 
range of parameters that needs to be scanned, at least for weaker couplings on smaller 
lattices. Obviously early stages of such work will be very much technical studies of 
the lattice theory. Nevertheless, we believe that the beginnings of first principles 
nonperturbative study of Af = 4 SYM are not so far off. As these progress, it will 
be interesting to compare our results to the on-going twisted supersymmetry lattice 
simulation studies that were initiated in [73]. 

5 Domain wall fermion lattice J\f = 1 SYM 

5.1 Domain wall fermions 

Lattice M = 1 super- Yang-Mills theory with GW fermions requires no fine-tuning. 
Domain wall fermions are a controllable approximation to GW fermions, and we have 
recently performed large scale simulations of the SU(2) theory [121-123]. We mea- 
sured the gaugino condensate, static potential, Creutz ratios and residual mass (a 
measure of explicit chiral symmetry breaking arising from the domain wall approxi- 
mation [124]). With this data we extrapolated the gaugino condensate to the chiral 
limit. We review some aspects of that study here. 

5.2 Introduction 

The only relevant or marginal operator allowed in a gauge invariant lattice formula- 
tion of pure M = 1 super- Yang-Mills [125] (SYM) with hypercubic symmetry is the 
gaugino mass term, as was emphasized long ago in the analysis of [56]. As above, 
GW lattice chiral symmetry protects against additive renormalizations of the gaugino 
mass in the continuum limit. Hence the desired continuum theory is obtained without 
fine-tuning of counterterms. 

The domain wall fermion (DWF) that we use originates from [97,126]. Properly 
speaking, it is GW only in the limit of infinite separation between the walls, L s — > oo. 
In extrapolations this is often traded for the residual mass m res , which is a measure 
of the explicit chiral symmetry breaking. 

Besides the absence of nonperturbative fine-tuning of the gaugino mass, DWF 
have the advantage that the fermion measure is positive and the square root of the 
determinant which enforces the Majorana condition is analytic with a phase that is 
independent of the gauge fields [58,60]. These three features are all lacking in the 
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Wilson fermion formulation that was applied in the only concerted lattice SYM effort 
to date, by the DESY-Miinster-Roma collaboration [61, 120, 127-130] and to a lesser 
extent Donini et al. [64,131,132]. (Recently, this program has been revived [133].) 
Our research is a continuation of the work of Fleming, Kogut and Vranas (FKV) [62] 
who first used DWF for studying Af = 1 SYM. Similar work has been initiated 
by Endres [134], with an extensive study appearing recently [135]. What sets the 
studies [121, 135] apart is that an extensive scan of the domain wall separation L s 
and measurement of the residual chiral symmetry breaking mass m res was done at 
different values of the bare lattice gauge coupling {(3 = 4/g 2 = 2.3 and 2.4 in our 
case) and spatial/temporal volumes (L 3 = 8 3 and 16 3 ; T = 16,32). This has allowed 
for chiral extrapolations (m rcs — > 0), and a preliminary view on what occurs as we 
take the continuum, theormodynamic limit (/3,L,T — > oo). 

It must be kept in mind that SYM does not have the Goldstone phenomena and 
the lightest states are the analogues of rj' and glueballs. The chiral symmetry that is 
broken in SYM is a discrete symmetry, Z 2 in the SU(2) study that we review here. 
The chiral regime is characterized by a gaugino mass that is small compared to the 
particle states of the theory, which should be of order l/r . Here r is the Sommer 
parameter [136], a measure of the dynamically generated length scale associated with 
confinement. Some of our results have rom res = 0.5, which is too large, but other 
results have r m rcs = 0.25, which ought to give chiral — hence supersymmetric — results 
to within 25%. 

5.3 Lattice formulation 

The lattice formulation that is used in this study has already been described by 
FKV [62,121]. It employs Shamir DWF [126] in the adjoint representation of SU(2) 
and the SU(2) fundamental plaquette Wilson gauge action. The Majorana condition 
is imposed through a square root on the fermion determinant, which as mentioned 
above is analytic and introduces no gauge field dependent sign ambiguity [58,60]. 

Lattice configurations were generated with a dynamical fermion mass rrif — 0, 
so that the finite size of the fifth dimension, parameterized by L s , was the sole in- 
frared regulator, through the corresponding additive mass correction m res , which is a 
measure of residual chiral symmetry breaking [124]. 

As was shown in the work by FKV, simulations performed at nonzero rrif, when 
extrapolated to mj = 0, give identical results to the rrif = simulations. Introducing 
rrif requires more simulations due to the rrif — > extrapolation that must be done. 



38 



Provided nif = simulations can be performed without crashing the fermion inverter, 
we believe this is preferable. 

Added confidence in the m,f — simulations comes from comparing our results to 
those of [135], as we discuss in the next section. 

5.4 Bare gaugino condensate 

A summary of all results obtained in [121] for the gaugino condensate (AA) is given 
in Tables [3] and H 

As mentioned above, confidence in the rrif = simulations comes from comparing 
our results to those of Endres [135]. He has performed simulations with mj ^ and 
extrapolated to m/ = 0, fitting to the linear function (AA)a 3 = cq + c\vrifa at fixed 
L s . His values of Cq for (3 = 2.3 on a 16 3 x 32 lattice compare well to the condensate 
we report in Table [3] at the L s values that can be matched, L s = 16, 24; cf. Table XII 
of [135] . Further comparisons will be discussed below; cf . Table [5J 

Measurements were conducted on large and small lattice volumes; it can be seen 
that in lattice units the finite-size dependence is mild or insignificant for (3 = 2.3 but 
quite noticeable for f3 = 2.4. This is sensible, given that (3 = 2.4 corresponds to a 
finer lattice spacing, and hence the physical volumes are smaller. Results below will 
show that the relative factor could be as large as 2 (cf. Table |6] vs. Table [7]). 

We measured the condensate at other values of L s using a sea-L s /valence-L s ap- 
proach. The condensate was measured using DWF with L™ 1 ' on top of dynamical 
lattices produced using a nearby L s s ea . Performing this for L s s ea values on either side 
of I/™ 1 - yields robust interpolated (I) results. 

We also used the results of the static potential study summarized in §5.5. II below 
to express m res and (AA) in terms of the Sommer scale r . Note that the (3 = 2.4 
value of m res r at L s = 48 indicates that the effective gaugino mass (which should 
be approximately equal to m res ) is roughly 1/4 the inverse Sommer scale, so that we 
are beginning to enter the chiral regime where supersymmetry is well approximated. 
On the other hand, it can be seen that m res r is unpleasantly large for (3 = 2.3 with 
L s < 32, and likewise the condensate in physical units is small compared to the 
(3 = 2.4 results. Clearly (3 = 2.3 is further away from the supersymmetric limit due 
to the coarser lattice. On the other hand it can be seen that the (3 = 2.4 data shows 
a marked volume dependence due to the smaller physical "box" that the states must 
squeeze into. 
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V x T 


L s 


m res a 


(AA)cr 


m res r 


(AA)r 3 


8 3 x 8 


16 


0.158(5) 


0.00711(7) 






8 3 x 32 


16 


0.181(3) 


0.00703(4) 


0.75(13) 


0.51(27) 


16 3 x 32 


16 


0.184(2) 


0.007051(5) 


0.668(10) 


0.337(11) 


8 3 x 32 


24 


0.1541(15) 


0.005112(8) 


0.610(97) 


0.32(15) 


16 x 32 


O A 

24 


r\ -1 p" f~* A ( "1 T\ 

0.1564(17) 


r\ aa r on 1 / r\\ 

0.005321(9) 


0.546(55) 


0.226(68) 


8 3 x 32 


32 


0.1319(12) 


0.004321(11) 


0.501(69) 


0.24(10) 


lb x 32 


32 


r\ -i H o /o\ 

0.143(2) 


r\ f\r\ A A C" /o\ 

0.00445(2) 


0.483(58) 


0.172(61) 


8 3 x 32 


40(1) 


0.1183(54) 


0.00383(3) 






8 3 x 16 


48 


0.1043(17) 


0.003563(20) 


0.361(31) 


0.148(37) 


8 3 x 32 


48 


0.1071(10) 


0.003551(11) 


0.409(31) 


0.198(45) 


8 3 x 32 


64 


0.08864(84) 


0.003164(10) 


0.300(35) 


0.122(42) 



Table 3: The gaugino condensate (AA) and residual mass m res for various lattice sizes 
and L s values, all at (3 = 2.3. The L s = 40 value, with an "(I)" after it, is obtained 
by the interpolation method described in the text. Values in units of the Sommer 
parameter ro are also shown, for those cases where the potential was measured (in 
particular, for all points that are included in the chiral extrapolation fit). The L s = 16 
data was not included in the linear chiral extrapolation fit, because these points had 
too much curvature (with respect to m rcs a) associated with them. 



V x T 


L s 




(AA)a 3 


m res r 


(AA)r 3 


8 3 x 32 


16 


0.080(2) 


0.004839(15) 


0.547(30) 


1.55(22) 


16 3 x 32 


16 


0.0969(8) 


0.00499(6) 


0.5355(66) 


0.842(25) 


8 3 x 32 


24 


0.0601(15) 


0.003293(17) 


0.417(26) 


1.10(18) 


16 3 x 32 


24 


0.0838(17) 


0.00389(8) 


0.385(35) 


0.38(10) 


16 3 x 32 


28(1) 


0.0721(33) 


0.003452(45) 






8 3 x 32 


32 


0.0486(12) 


0.00269(2) 


0.296(15) 


0.61(08) 


16 3 x 32 


32 


0.0653(15) 


0.003330(12) 


0.313(33) 


0.37(11) 


8 3 x 32 


40(1) 


0.0390(24) 


0.00234(8) 






8 3 x 32 


48 


0.0328(9) 


0.002165(18) 


0.224(17) 


0.69(15) 



Table 4: Results similar to Table [HI except that these are for (3 = 2.4. 
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Endres/Extrap. 


Us (8 3 x 16) 


Us (16 3 x 32) 




(AA)a 3 (16 3 x 32) 


(AA)a 3 


(AA)a 3 


16 


0.0070544(51) 


0.00703(4) 


0.007051(5) 


20 


0.0058979(55) 


— 


— 


24 


0.0051697(49) 


0.005112(8) 


0.005321(9) 


28 


0.0046770(51) 






32 


0.00432(2) 


0.004321(11) 


0.00445(2) 


40 


0.00381(2) 


0.00383(3) 




48 


0.00346(3) 


0.003551(11) 





Table 5: Column 2 contains extrapolation results of Endres for L s = 16, 20, 24, 28 
with (3 = 2.3 on the 16 3 x 32 lattice. Also in this column are the results we obtain for 
L s = 32, 40, 48 from his data. Here, we model the behavior of c\ from (15. ip as given 
by (15.21) . extracted from his results. We then use his rrif = 0.02 results at these L s 
values to estimate the rrif = condensate. In columns 3 and 4 we give our simulation 
results for comparison. 

5.4.1 Combined results for (3 = 2.3 

Here we compar in detail our results to those of Endres [135]. At each L s he has fit 
rrif results to the form: 

(AA)a 3 = Co + ci rrifa. (5.1) 

What is interesting about the results of his Table XI is that the coefficient C\ shows 
a regular pattern with respect to L s : 

c x {L s ) = 0.0828(1) + 0.00025(7) x (L s - 28). (5.2) 

Here, the numbers in parentheses represent our estimate of the error in this formula, 
based on his results. From this formula we can use his measurements of (AA)a 3 
at rrif = 0.02 for L s = 32,40,48 to obtain predictions of the rrif = condensate. 
These are shown, combined with his extrapolation values at L s = 16, 20, 24, 28 in 
Table [5j Comparing to our Table [3J results of which we reproduce in the two right- 
most columns of Table (H we find that there is reasonable consistency. We will use 
these combined results in the chiral extrapolation fit below. 
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5.5 Gluonic observables 



5.5.1 Static potential 

The static potential was obtained by measuring Wilson loops with one side of length 
t in the temporal direction, according to standard methods. Having obtained V{r)a 
from fitting the exponential decay in time, we fit the data to the standard form 



This approach to the determination of r$/a has some sensitivity to the range of radii 
that is fit, and obviously depends on what form we assume for V(r)a. 

The results of our static potential fit are presented in Tables |6] and [71 For the 
L = 16 results, the fits were also done using the same set of Wilson loops as in the 
L = 8 case, denoted "L = 8 method" , so that dependence choice of Wilson loops 
could be controlled for, and therefore ruled out as a spurious source of finite size 
dependence. 

These results are to be compared with Table 1 of [134] or Table IX of [135]. 
There, a nonzero fermion mass was used different fit ranges for r, t were employed. 
In particular, our fits include the points r/a = 1, V2, which have very small errors 
and thus strongly influence the fit. Thus for L s = 16 we have also performed a fit 
with r m j n = v3 as was done in [134], as can be seen in the second L s = 16 entry for 
(3 = 2.3, 2.4 in Tables [6] and [3 Our results for the fit quantities at (3 = 2.3 are in good 
agreement once this restriction is imposed. For our other L s values we have far fewer 
samples, as larger L s comes at greater computing cost. The degradation of statistical 
errors that results if we exclude the r = 1, \/2 points is unacceptable, which is why 
we do not quote results with the same r min = \/3 as [134] for the other L s values. On 
the other hand, it can be seen from the L s = 16 results that the choice of r m ; n only 
has a 10% effect on the r /a estimate, so that the choice of r min is not crucial to the 
broad picture that we are after. The (3 = 2.4 results are also in reasonable agreement 
with [134], comparing to the numbers we obtain at 16 3 x 32, L s = 16, with r min = \^3, 
and keeping in mind the nonzero mj in [134]. 

Above, we have used the results of Tables [6] and [7J to scale the residual mass and 
condensate to r units. (Note that the r /a values with identical lattice parameters 



V{r)a = V a + aa 2 (r/a) — . 

r/a 

We obtain the Sommer parameter r Q /a from this fit, using the formula 



(5.3) 




(5.4) 
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8 3 x 32 


32 


0.748(82) 


0.087(27) 


0.400(55) 


3.80(52) 


1.250(55) 


L = 8 


16 3 x 32 


32 


0.745(90) 


0.109(31) 


0.412(59) 


3.38(40) 


1.238(59) 


L = 16 


16 3 x 32 


32 


0.635(50) 


0.146(17) 


0.338(33) 


3.00(14) 


1.312(33) 


L = 8 


8 3 x 16 


48 


0.706(68) 


0.107(22) 


0.372(47) 


3.46(29) 


1.278(47) 


L = 8 


8 3 x 32 


48 


0.768(47) 


0.085(15) 


0.414(33) 


3.82(29) 


1.236(33) 


L = 8 


8 3 x 32 


64 


0.680(94) 


0.113(32) 


0.353(63) 


3.38(39) 


1.297(63) 


L = 8 



Table 6: Gluonic observables obtained from the static potential for (3 = 2.3. The 
quantity ro is physical, so that the ratio r^/a provides a measure of the lattice spacing. 



V x T 


L. 


V a 


2 

era 


a 


ro/a 


arl 


method 


8 3 x 32 


16 


0.617(11) 


0.0292(30) 


0.2857(83) 


6.84(33) 


1.3643(83) 


L = 8 


16 3 x 32 


16 


0.5846(32) 


0.04531(91) 


0.2659(24) 


5.526(51) 


1.3841(24) 


L= 16 


16 3 x 32 


16 


0.537(11) 


0.0554(20) 


0.219(15) 


5.083(67) 


1.431(15) 


L = 16, r > V% 


8 3 x 32 


24 


0.636(12) 


0.0280(33) 


0.2997(88) 


6.94(39) 


1.3503(88) 


L = 8 


16 3 x 32 


24 


0.579(40) 


0.065(13) 


0.272(27) 


4.60(41) 


1.378(27) 


L = 16 


8 3 x 32 


32 


0.609(12) 


0.0369(36) 


0.2809(90) 


6.09(28) 


1.369(90) 


L = 8 


16 3 x 32 


32 


0.611(43) 


0.059(13) 


0.295(29) 


4.79(50) 


1.355(29) 


L= 16 


8 3 x 32 


48 


0.648(15) 


0.0288(44) 


0.309(11) 


6.83(49) 


1.341(11) 


L = 8 



Table 7: Gluonic observables obtained from the static potential for (3 = 2.4. We note 
that r /a shows significant volume dependence for L s = 32. Comparing to Table [61 
the L, L s —>■ oo trend seems to indicate a lattice spacing that is slightly more 1/2 of 
the 13 = 2.3 value, in agreement with what we estimate using two-loop renormalization 
of the gauge coupling. 
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were used in this procedure, rather than a uniform r /a value across all m res and (AA).) 
With the string tension in hand, we now see that the energy scale of confinement 
a/ otq « 1.4 lies above the explicit chiral symmetry breaking scale m res r by a factor 
of 1.8 to 3.8 for (3 = 2.3, and 3.0 to 5.2 for (3 = 2.4. This is consistent with the 
observation that the string tension results in Tables [6] and [7| are insensitive to the 
range of L s values displayed there, when expressed in physical units (otq). That 
is, confinement dynamics are to a good approximation decoupled from the explicit 
chiral symmetry breaking. Since the lowest lying excitations of SYM are glueballs 
and superpartners, the gap associated with confinement should also decouple these 
states from the explicit chiral symmetry breaking. Thus it appears that we in the 
regime where the spectrum reflects supersymmetry, modulo lattice artifacts, and it 
will be quite interesting to examine the spectrum in order to check whether or not 
this is true — something we will do in future work. 

At this point we make some further remarks regarding finite size effects, since 
the volumes are for several points small. Because there are no Goldstone modes the 
situation is very different from QCD, where one has to watch out for the pion Compton 
wavelength. Instead, here we have to watch out for the rf Compton wavelength, which 
will be of order r . 

For the 16 3 x 32 lattices at (3 — 2.3, we have r^/a = 3.0 to 3.5, which gives a lattice 
extension 4.6 to 5.3 tq. If the glueballs and rj' have wavelength tq, we should be 
safe from finite volume effects. On the 8 3 lattices there is a bigger problem, but still 
the lattice extent is always at least 2r for the points that are used in our chiral 
extrapolation for (3 = 2.3. In fact, our results in Table II show only a mild volume 
dependence for bare lattice quantities, confirming the arguments we have just made. 

The volume dependence for (3 = 2.4 is more significant, as we have discussed 
above. For the 16 3 x 32 lattices we have lattice extent 3.3r to 3.5r in space-like 
directions. For the 8 3 lattices we have lattice extent 1.2ro, which is too small. This 
likely explains why the L = 8 versusy L = 16 results are so different for the measured 
quantites. In fact, it may only be the L = 8 results that are significantly impacted 
by finite size effects; the L = 16 lattices, having 3.3r to 3.5r , are perhaps large 
enough. Evidence in that direction is given by the insensitivity of the (3 = 2.3 results 
for lattices of this size, in units of tq. 

In the case of (3 = 2.4, further simulations should eventually be done on lat- 
tices larger than 16 3 x 32, but it is a major undertaking that is beyond our present 
capabilities, which are at present state of the art. 
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5.6 Extrapolation of the gaugino condensate 

In our previous work, we presented only a linear extrapolation of the condensate in 
m res , to be discussed in §5.6.11 Here we will also present some new results where we 
perform nonlinear chiral extrapolations based on the functional forms that Endres 
has assumed in [135]. It will be seen that although our simulations are in agreement, 
different assumptions about the functional form of the chiral extrapolation lead to 
m res — > estimates that differ by factors as large as 3, as has already been noted 
in [135]. 

5.6.1 Linear extrapolation 

One important question is the size of L s necessary to get into the linear regime where 

(AA)a 3 = c + C]m res a (5.5) 

is a good approximation. Obviously, this serves as an indicator of where we need to be 
in order to have SYM well-approximated. Thus, the measurement of (AA) vs. m res is 
an important benchmark for determining the regime in which other SYM phenomena 
can be studied with the DWF lattice approach. Another question is the extent to 
which cq,! are sensitive to finite spacetime volume (V4 = VxT = L 3 xT in our 
notation). Interestingly, we find that most of the volume dependence is absorbed into 
be seen from the dashed lines in Figs. [2] and [31 To a good approximation 
the 8 3 x 32 and 16 3 x 32 lattice data lie on the same line. The smaller value of m res 
on the smaller lattice is most likely due to a smaller density of near-zero modes. The 
linear chiral extrapolation (m rcs — > 0) of (AA)a 3 obtained from the fit to (15. 5p is given 
in Table [H A feel for the sensitivity to the fitted range of L s can be seen from the two 
results we provide for (3 = 2.3, which differ by the minimum L s that was included. In 
fact, the quality of the L s > 16 fit is very poor due to nonlinear dependence on m res 
that enters at L s = 24, as can also be seen from Fig. [21 

5.6.2 Nonlinear extrapolation 

We now follow Endres [135] and fit the data to the form 

(AA)a 3 = b + exp(-6 2 L s ). (5.6) 

All (3 = 2.3 results of Tables [3] and [5] except our L s = 16 results for 8 3 x 8 and 8 3 x 32 
lattices are included. We then exclude L s = 16 data, and then L s = 16, 20 data. The 
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(3 


L s range 


X 2 /d.o.f. 


Co 




2.3 


24-64 


136 


0.00026(25) 


0.0316(20) 


2.3 


32-64 


30.0 


0.00086(17) 


0.0258(15) 


2.4 


24-48 


19.5 


0.00098(13) 


0.0364(23) 



Table 8: Fit results for the linear extrapolation ( 15. 51) of the gaugino condensate in 
m res , depending upon the range of L s values used. For (3 = 2.3, the quality of the 
L s > 16 fit is very poor due to nonlinear dependence on m res that enters at L s = 24, 
as can also be seen from Fig. [2j 




Figure 2: Condensate vs. m res for (3 = 2.3, in bare lattice units. Dashed lines show the 
two linear fits (differing by the minimum L s included) . Extrapolated values together 
with fit errors are shown at m rcs = 0. The Solid line combines results from fitting 
functions (15. 6p and (15.71) . 



46 




Figure 3: Condensate vs. m res for (3 = 2.4, in bare lattice units. The dashed line shows 
the linear fit. The extrapolated value together with fit error is shown at m res = 0. 
Solid lines combine results from fitting functions ( 15. 6ft and (15.7ft . 



results of the fit are presented in Table [91 The chiral limit of (AA)a 3 is given by b^. 
It can be seen that (AA)a 3 = 0.0027(2) is a good characterization of the result. We 
note that this is much larger than the value obtained by linear extrapolation in m res . 

In Table fTUl we fit our (3 = 2.4 to the same formula. Due to the large finite size 
effects at this finer lattice spacing, we have separately fit the L = 8 and L = 16 data. 
The dependence on L is also strong in the chiral extrapolation of the condensate, bo. 
We would need an additional L value, such as L = 12, in order to obtain an estimate 
of the condensate in the L — > oo limit. Nevertheless we see that the trend is for the 
condensate to increase with L, and we obtain an estimate of (AA)a 3 ^ 0.0029 in the 
chiral limit. We argued above that the L = 16 lattices, having 3.3r to 3.5r , are 
perhaps large enough to give a good approximation of the L — > oo behavior. If that 
is true, then we obtain an estimate for the (3 = 2.4 chiral limit of (AA)a 3 = 0.0029(4). 
Note that this is significantly larger than the result of linear extrapolation in m res . 

Endres has argued that the residual mass should be fit to 

m res a = — (a exp(-aiL s ) + a 2 ), (5.7) 
a functional form that was derived by transfer matrix methods in Eq. (3.7) of [137]. 
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L s range 


X 2 /d.o.f. 


bo 


h 


b 2 


16-64 


24.5 


0.002796(96) 


0.0950(24) 


0.0209(27) 


20-64 


26.6 


0.00270(16) 


0.0901(53) 


0.0171(51) 


24-64 


28.6 


0.00275(20) 


0.093(10) 


0.0193(78) 



Table 9: Fit of equation ( 15. 6ft for (3 = 2.3. 



L 


X 2 /d.o.f. 


bo 




b 2 


8 


0.699 


0.001848(71) 


0.0792(82) 


0.0342(62) 


16 


3.27 


0.00292(40) 


0.084(51) 


0.058(49) 



Table 10: Fit of equation (15.61) for (3 = 2.4, using all data from Table H] for each 
L = 8,16. 

Endres has shown that (15. 7p is consistent with his rrif — 0.02 data. Here we have 
performed a fit to this function using our m/ = data. The results are given in 
Table [TTJ In order to obtain stable fits, it was necessary to separate the m res data for 
(3 = 2.4 according the two L values. One can see that the finite size dependence of 
the fit parameters is significant in that case. 

Combining the results of the fit functions (15.61) and (15.71) . we obtain a map L s — > 
(m res a, (AA)a 3 ) for arbitrary L s . With this we obtain the nonlinear fit curves that 
appear in Figs. [2] and [3] as solid lines. 

5.6.3 Discussion 

The message is that one needs data at smaller m res . The nonlinear fits are motivated 
by known features of the DWF discretization and are more likely to be trustworthy 
at the large m ies values that we are working with. Still, it would be greatly reassuring 



p 


L 


L s range 


X 2 /d.o.f. 


a 


ai 


a-2 


2.3 


all 


16-64 


4.47 


-5.71(22) 


0.0308(58) 


6.44(38) 


2.4 


8 


16-48 


0.250 


-1.64(56) 


0.102(24) 


1.596(28) 


2.4 


16 


16-32 


0.237 


-20(17) 


0.226(55) 


2.097(34) 



Table 11: Fit of m res data to ( 15.71) . 
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to have simulations at small enough m res that the linear extrapolation agrees with 
these formulae. At a minimum, simulations at larger L s should be performed as a 
test of the predictions made by the solid curves in Figs. [2] and [31 

According to Table [TJ the value of the lattice spacing a is smaller on the j3 = 
2.4 lattice. Thus it is surprising that the extrapolated values of (AA)a 3 is as large 
for j3 = 2.4 as it is for (3 = 2.3. A plausible interpretion is that there are larger 
renormalizations of the condensate on a finer lattice (J3 = 2.4), a hypothesis that we 
are preparing to test with nonperturbative renormalization [138, 139] in an upcoming 
study. 

5.7 Outlook 

Future work that is envisioned aims to develop a deeper understanding of the con- 
figurations that are responsible for generating the nonzero gaugino condensate. In 
particular, we would like to elucidate the continuum picture on the cylinder R 3 x S 1 , 
where it is monopoles and "KK monopoles" that combine to yield the infinite volume 
value [140] . It is already known from spectral flow techniques that fractional topolog- 
ical charge plays an important role in the zeromodes of adjoint lattice fermions [141], 
and we would like to explicitly connect that with the semiclassical configurations 
that have been suggested by continuum methods. It would also be interesting to 
study questions of center dominance and abelian dominance in LSYM, following the 
line of research that has been pursued in the "pure glue" (i.e. quenched) SU(2) the- 
ory [142-145]. 

Note that rather large L s values were required in order to get m res ro ~ 1/4. To 
improve the situation we envision switching to simulations with modified versions of 
DWF that have superior chiral behavior [146,147]. This should clear up uncertainties 
in the chiral extrapolation that are apparent from the results in §5.61 above. 

6 Conclusions 

As algorithms and hardware have advanced, so have the prospects for simulating 
four- dimensional supersymmetric gauge theories. We have illustrated in three differ- 
ent types of models how GW chiral symmetry and exact lattice flavor symmetries can 
be used to limit the counterterms that must be adjusted to get the desired continuum 
limit. Furthermore, it was seen that in most cases only bosonic counterterms need 
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to be adjusted. Multicanonical reweighting techniques provide a viable method for 
accomplishing this on small lattices in the near term. On the other hand, the Af = 1 
SYM study that we have concluded with illustrates the magnitude of the challenge in 
simulating larger lattices. In fact, that study consumed 30 million IBM BlueGene/L 
node hours. Nevertheless, continuing progress in the program of simulating super- 
symmetric gauge theories is to be expected and we look forward to reporting our 
advances in this field as further results, especially in SQCD and Af = 4 SYM, are 
obtained. 
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Appendices 

A Super-QCD continuum theory 

Here we describe the supersymmetric continuum theory with massless "quarks" and 
scalar partners, "squarks." Since lattice studies of SQCD are intended for under- 
standing new strong interactions that are associated with spontaneous supersymme- 
try breaking and its transmission to the visible sector, these terms are only used by 
way of analogy; the SU(N) gauge symmetry is in addition to the Standard Model 
gauge group, and is strongly coupled at scales of a TeV or perhaps much more. 
The Minkowski space Lagrangian is 

C = -\ Tr G^G^ - 2i Tr Xa^V^X + Tr D 2 - (V^Wp 

-iXp^V^Xp + F ] p F p + iV2gp^Xx P - iV2gx p Xp 

+gp ] Dp+ (p -> q). (A.l) 

Here we mostly follow the two-component conventions of [91], except that we use /i, v 
for spacetime indices and write the "gluon" field strength G^ u . We have auxiliary 
scalar fields D = D a t a and F p which allow for supersymmetry without imposing the 
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equations of motion (off-shell supersymmetry). They are easily integrated out using 
their equations of motion. Flavor indices have been suppressed, and our convention for 
the generators is Tr t a t b = (l/2)S ab . As in £J21 we use a matrix notation in SU(N) color 
space, so that Xp is supposed to represent an iV-component row vector. Covariant 
derivatives X> M take the usual form. The Yukawa terms can be related to those of the 
main text by the redefinition 

A -f iX, X — > -iX. (A.2) 

Once this is done, it is clear that the CP and S symmetries that have been im- 
posed in §|2]are also symmetries of the continuum theory (lA.lj) . The supersymmetry 
transformations under which (lA.ip is invariant are, for the matter fields, 



5F p = iV2^V^ + i2gfXp, (p-g), (A.3) 

and for the gauge fields, 

5^ = -iXa^ + i^X, 5^X = a^^G^ + %D, 

5^D = -i^V^X - (A.4) 

where £ is a Grassmann transformation parameter. It is straightforward to continue 
to imaginary time and the Euclidean version of the action that we use in £j2j 



References 

[1] J. Giedt, "Deconstruction and other approaches to supersymmetric lattice field 
theories," Int. J. Mod. Phys. A21 (2006) 3039, jarXiv : hep-lat/0602007] . 

[2] E. Witten, 'An SU(2) anomaly," Phys. Lett. B117 (1982) 324. 

[3] M. Hayakawa and H. Suzuki, "Gauge anomaly associated to the Majorana 
fermion in 8k + 1 dimensions," Prog. Theor. Phys. 115 (2006) 1129, 
(arXiv : hep-th/0601026| . 

[4] A. Casher and Y. Shamir, "Feynman rules for non-perturbative sectors and 
anomalous supersymmetry ward identities." arXiv:hep-th/9908074. 



51 



[5] N. Seiberg, "The Power of holomorphy: Exact results in 4-D SUSY field 
theories," arxiv:hep-th/9408013. Talk given at Particles, Strings, and 
Cosmology (PASCOS 1994), Syracuse, New York, 19-24 May, 1994. 



[6] V. S. Kaplunovsky and J. Louis, "Model independent analysis of soft terms in 
effective supergravity and in string theory," Phys. Lett. B306 (1993) 269, 
|arXiv : hep-th/9303040| . 

[7] A. Brignole, L. E. Ibanez and C. Munoz, "Towards a theory of soft terms for 
the supersymmetric standard model," Nucl. Phys. B422 (1994) 125, 
|arXiv : hep-ph/930827l] . [Erratum-ibid. B 436 (1995) 747]. 

[8] A. G. Cohen, T. S. Roy, and M. Schmaltz, "Hidden sector renormalization of 
MSSM scalar masses," JEEP 02 (2007) 027, |hep-ph/06 12100] . 

[9] T. S. Roy and M. Schmaltz, "A hidden solution to the fi/B^ problem in gauge 
mediation," Phys. Rev. D77 (2008) 095008, jarXiv : 07087359 3] . 

[10] H. P. Nilles, "Gaugino condensation and SUSY breakdown," |hep-th/0402022[ 

[11] M. K. Gaillard and B. D. Nelson, "Kahler stabilized, modular invariant 
heterotic string models," Int. J. Mod. Phys. A22 (2007) 1451, 
|hep-th/0703227| . 

[12] G. F. Giudice and R. Rattazzi, "Theories with gauge-mediated 

supersymmetry breaking," Phys. Rept. 322 (1999) 419, [hep-ph/9801271| . 

[13] A. G. Cohen, D. B. Kaplan, and A. E. Nelson, "The more minimal 
supersymmetric standard model," Phys. Lett. B388 (1996) 588-598, 
|hep-ph/9607394| . 

[14] N. Arkani-Hamed, M. A. Luty, and J. Terning, "Composite quarks and 

leptons from dynamical supersymmetry breaking without messengers," Phys. 
Rev. D58 (1998) 015004, hep-ph/97 123891 . 

[15] M. A. Luty and J. Terning, "Improved single sector supersymmetry breaking," 
Phys. Rev. D62 (2000) 075006, hep-ph/98 12290] . 

[16] R. R. Volkas and G. C. Joshi, "Supersymmetric composite models," Phys. 
Rept. 159 (1988) 303. 



52 



[17] M. J. Strassler, "Generating a fermion mass hierarchy in a composite 
supersymmetric standard model," Phys. Lett. B376 (1996) 119, 
|hep-ph/9510342| . 

[18] A. E. Nelson and M. J. Strassler, "A realistic supersymmetric model with 
composite quarks," Phys. Rev. D56 (1997) 4226, |arXiv : hep-ph/9607362] . 

[19] M. Gabella, T. Gherghetta, and J. Giedt, 'A gravity dual and LHC study of 
single-sector supersymmetry breaking," Phys. Rev. D76 (2007) 055001, 
jarXiv: 0704. 3571] . 

[20] T. Gherghetta and J. Giedt, "Bulk fields in AdS(5) from probe D7 branes," 
Phys. Rev. D74 (2006) 066007, hep-th/0605212] . 

[21] J. Giedt, "Completion of standard model-like embeddings," Ann. Phys. 289 

(2001) 251, |hep-th/0009"T04] . 

[22] J. Giedt, "Spectra in standard-like Z(3) orbifold models," Annals Phys. 297 

(2002) 67-126, |hep-th/0108244] . 

[23] J. Giedt, "Z(3) orbifolds of the SO(32) heterotic string: 1 Wilson line 
embeddings," Nucl. Phys. B671 (2003) 133-147, |hep-th/030l"232] . 

[24] J. Giedt, "Lattice supersymmetry and string phenomenology," in proceedings 
of 2nd International Conference on String Phenomenology 2003, Durham, 
England, 29 Jul - 4 Aug 2003 (2003) |hep-ph/0309344] . 

[25] J. Giedt, "The KM phase in semi-realistic heterotic orbifold models," Nucl. 
Phys. B595 (2001) 3-32, | hep-ph/0007193j . 

[26] M. K. Gaillard and J. Giedt, "A D-moduli problem?," Phys. Lett. B479 
(2000) 308-314, hep -ph/00012191 . 

[27] J. Giedt, "CP violation and moduli stabilization in heterotic models," Mod. 
Phys. Lett. A17 (2002) 1465-1473, |hep-ph/0204017] . 

[28] J. Giedt, "Optical unification," Mod. Phys. Lett. A18 (2003) 1625-1633, 
|hep-ph /0205224| . 



53 



[29] J. Giedt, "D-moduli stabilization," in proceedings 1st International 

Conference on String Phenomenology (SP2002), Oxford, England, 6-11 Jul 
2002 (2002) [he p-ph/0208004] . 

[30] M. K. Gaillard, J. Giedt, and A. L. Mints, "Modular invariant gaugino 

condensation in the presence of an anomalous U(l)," Nucl. Phys. B700 (2004) 
205-270, |ne p-th/0312125| . 

[31] J. Giedt, "Lack of trinification in Z(3) orbifolds of the SO(32) heterotic 
string," Mod. Phys. Lett. A20 (2005) 2369-2375, |hep-ph/040220T1 . 

[32] J. Giedt, G. L. Kane, P. Langacker, and B. D. Nelson, "Massive neutrinos and 
(heterotic) string theory," Phys. Rev. D71 (2005) 115013, |hep-th/0502032] . 

[33] J. R. Ellis, J. Giedt, O. Lebedev, K. Olive, and M. Srednicki, "Against 
Tachyophobia," Phys. Rev. D78 (2008) 075006, [arXiv: 0806. 3648] . 

[34] T. Bhattacharya, M. R. Martin, and E. Poppitz, "Chiral lattice gauge theories 
from warped domain walls and Ginsparg- Wilson fermions," Phys. Rev. D74 



(2006) 085028, |hep-lat/0605003|. 

[35] J. Giedt and E. Poppitz, "Chiral Lattice Gauge Theories and The Strong 
Coupling Dynamics of a Yukawa-Higgs Model with Ginsparg- Wilson 



Fermions," JEEP 10 (2007) 076, |hep -lat/0701 004|. 

[36] E. Poppitz and Y. Shang, "Lattice chirality and the decoupling of mirror 
fermions," JEEP 08 (2007) 081, jarXiv : 0706TT04 3] . 

[37] E. Poppitz and Y. Shang, '"Light from chaos' in two dimensions," Int. J. 
Mod. Phys. A23 (2008) 4545-4556, jarXiv: 0801 . 0587] . 

[38] E. Poppitz and Y. Shang, "Lattice chirality, anomaly matching, and more on 
the (non) decoupling of mirror fermions," arXiv: 0901 .3402, 

[39] S. Catterall and T. Wiseman, "Black hole thermodynamics from simulations 
of lattice Yang-Mills theory," Phys. Rev. D78 (2008) 041502, 
jarXiv: 0803. 4273] . 

[40] S. Catterall and T. Wiseman, "Gauge-gravity duality — super Yang-Mills 
quantum mechanics," PoS LAT2007 (2007) 051, [ arXiv: 0709 . 3497] . 



54 



[41] S. Catterall and T. Wiseman, "Towards lattice simulation of the gauge theory 
duals to black holes and hot strings," JEEP 12 (2007) 104, 
jarXiv: 0706. 3518] . 

[42] M. Hanada, J. Nishimura, and S. Takeuchi, "Non-lattice simulation for 

supersymmetric gauge theories in one dimension," Phys. Rev. Lett. 99 (2007) 
161602, |arXiv:0706.1647j . 

[43] T. Banks, W. Fischler, S. H. Shenker and L. Susskind, "M theory as a matrix 
model: A conjecture," Phys. Rev. D55 (1997) 5112, |arXiv:hep-th/96 10043] . 

[44] N. Ishibashi, H. Kawai, Y. Kitazawa and A. Tsuchiya, "A large-N reduced 
model as superstring," Nucl. Phys. B498 (1997) 467, 
|arXiv : hep-th/96 12115] . 

[45] J. M. Maldacena, "The large N limit of superconformal field theories and 
supergravity," Adv. Theor. Math. Phys. 2 (1998) 231, 
|arXiv:hep-th/97lT200] . 

[46] S. S. Gubser, I. R. Klebanov and A. M. Polyakov, "Gauge theory correlators 
from non-critical string theory," Phys. Lett. B428 (1998) 105, 
|arXiv:hep-th/9802109| . 

[47] E. Witten, "Anti-de Sitter space and holography," Adv. Theor. Math. Phys. 2 
(1998) 253, |arXiv : hep-th/9802150l . 

[48] R. R. Metsaev, "Type IIB Green-Schwarz superstring in plane wave 
Ramond-Ramond background," Nucl. Phys. B625 (2002) 70, 
|arXiv : hep-th/01 12044| . 

[49] R. R. Metsaev and A. A. Tseytlin, "Exactly solvable model of superstring in 
plane wave Ramond-Ramond background," Phys. Rev. D65 (2002) 126004, 
[arXiv : hep-th/0202109] . 

[50] M. Blau, J. Figueroa-O'Farrill, C. Hull and G. Papadopoulos, "A new 
maximally supersymmetric background of IIB superstring theory," JEEP 
0201 (2002) 047, |arXiv: hep-th/01 10242 . 



55 



[51] D. Berenstein, J. M. Maldacena and H. Nastase, "Strings in flat space and PP 
waves from Af = 4 super- Yang-Mills," JEEP 0204 (2002) 013, 
[arXiv : hep-th/0202021 . 

[52] P. H. Dondi and H. Nicolai, "Lattice supersymmetry," Nuovo Cim. A41 
(1977) 1. 

[53] K. Fujikawa, "Supersymmetry on the lattice and the Leibniz rule," Nucl. 
Phys. B636 (2002) 80-98, lhep-th/0205095] . 

[54] J. Giedt and E. Poppitz, "Lattice supersymmetry, superfields and 
renormalization," JEEP 0409 (2004) 029, |hep-th/0407135| . 

[55] P. H. Ginsparg and K. G. Wilson Phys. Rev. D25 (1982) 2649. 

[56] G. Curci and G. Veneziano, "Supersymmetry and the lattice: A 
reconciliation?," Nucl. Phys. B292 (1987) 555. 

[57] N. Maru and J. Nishimura, "Lattice formulation of supersymmetric 

Yang-Mills theories without fine-tuning," Int. J. Mod. Phys. A13 (1998) 2841, 
[arXiv : hep-th/9705152] . 

[58] H. Neuberger, "Vector like gauge theories with almost massless fermions on 
the lattice," Phys. Rev. D57 (1998) 5417, [ hep-lat/9710089 . 

[59] J. Nishimura, "Four-dimensional M = 1 supersymmetric Yang-Mills theory on 
the lattice without fine-tuning," Phys. Lett. B406 (1997) 215-218, 
|arXiv:hep-lat/9701013| . 

[60] D. B. Kaplan and M. Schmaltz, "Supersymmetric Yang-Mills theories from 
domain wall fermions," Chm. J. Phys. 38 (2000) 543, 
|arXiv : hep-lat/0002030| . 

[61] I. Campos, R. Kirchner, I. Montvay, J. Westphalen, A. Feo, S. Luckmann, 
G. Minister, K. Spanderan [DESY-Munster Collaboration], "Monte Carlo 
simulation of SU(2) Yang-Mills theory with light gluinos," Eur. Phys. J. Cll 
(1999) 507, |arXiv : hep-lat/9903014] . 

[62] G. T. Fleming, J. B. Kogut and P. M. Vranas, "Super- Yang-Mills on the 
lattice with domain wall fermions," Phys. Rev. D64 (2001) 034510, 
|arXiv : hep-lat /00 08009 1 . 



56 



[63] DESY-Munster Collaboration, R. Kirchner, S. Luckmann, I. Montvay, 
K. Spanderen, and J. Westphalen, "Numerical simulation of dynamical 
gluinos: Experience with a multi-bosonic algorithm and first results," Nucl. 
Phys. Proc. Suppl. 73 (1999) 828, [ jarXiv : hep-lat/9808024] . 

[64] A. Donini, M. Guagnelli, P. Hernandez, and A. Vladikas, "Towards M = 1 
super- Yang-Mills on the lattice," Nucl. Phys. B523 (1998) 529, 
|hep-lat/9710065| . 

[65] D. B. Kaplan, E. Katz and M. Unsal, "Supersymmetry on a spatial lattice," 
JEEP 0305 (2003) 037, h ep-lat/0206019| . 

[66] A. G. Cohen, D. B. Kaplan, E. Katz and M. Unsal, "Supersymmetry on a 
euclidean spacetime lattice. II: Target theories with eight supercharges," 
JEEP 0312 (2003) 031, |hep-lat/0307012| . 

[67] S. Catterall, "Lattice formulation of M = 4 super-yang-mills theory," JEEP 
0506 (2005) 027, [ arXiv : hep-lat/050303"6 . 

[68] D. B. Kaplan and M. Unsal, "A Euclidean lattice construction of 

supersymmetric Yang-Mills theories with sixteen supercharges," JEEP 09 
(2005) 042, |arXiv : hep-lat/0503039l . 

[69] S. Catterall, "From Twisted Supersymmetry to Orbifold Lattices," JEEP 01 
(2008) 048, jarXiv: 0712. 2532] . 

[70] J. Giedt, "Non-positive fermion determinants in lattice supersymmetry," Nucl. 
Phys. B668 (2003) 138, |hep-lat/0304006| . 

[71] J. Giedt, "The fermion determinant in (4,4) 2d lattice super- Yang-Mills," 
Nucl. Phys. B674 (2003) 259, |hep-lat/0307024| . 

[72] J. Giedt, "Deconstruction, 2d lattice super- Yang-Mills, and the dynamical 
lattice spacing," |hep-lat/0405021 

[73] S. Catterall, "First results from simulations of supersymmetric lattices," 
JEEP 01 (2009) 040, [arXiv: 081 17120 3]. 

[74] I. Kanamori, "Vacuum energy of two-dimensional N = (2, 2) super Yang-Mills 
theory," larXiv: 0902 .28761 



57 



[75] I. Kanamori and H. Suzuki, "Some physics of the two-dimensional M = (2, 2) 
supersymmetric Yang-Mills theory: Lattice Monte Carlo study," Phys. Lett. 
B672 (2009) 307, jarXiv: 0811 . 2851] . 

[76] S. Catterall, "Dirac-Kahler fermions and exact lattice supersymmetry," PoS 
LAT2005 (2006) 006, larXiv:hep-lat/0509136] . 

[77] J. Giedt, "Advances and applications of lattice supersymmetry," PoS 



LAT2006 (2006) 008, | arXiv : hep -lat/0701 006 1 . 

J. W. Elliott, J. Giedt and G. D. Moore, "Lattice four-dimensional N = 4 
SYM is practical," Phys. Rev. D78 (2008) 081701, j arXiv: 0806. 0013] . 

[79] M. G. Endres and D. B. Kaplan, "Lattice formulation of (2,2) supersymmetric 
gauge theories with matter fields," JEEP 10 (2006) 076, |hep-lat/0604012| . 

[80] J. Giedt, "Quiver lattice supersymmetric matter, D1/D5 branes and 
AdS(3)/CFT(2)," |hep-lat/0605004| 

M. Falcioni, E. Marinari, M. L. Paciello, G. Parisi and B. Taglienti, "Complex 
zeros in the partition function of the four-dimensional SU(2) lattice gauge 
model," Phys. Lett. B108 (1982) 331. 

[82] A. M. Ferrenberg and R. H. Swendsen, "New Monte Carlo technique for 
studying phase transitions," Phys. Rev. Lett. 61 (1988) 2635. 

[83] A. M. Ferrenberg and R. H. Swendsen, "Optimized Monte Carlo analysis," 
Phys. Rev. Lett. 63 (1989) 1195. 

[84] K. Binder, "Finite size scaling analysis of Ising model block distribution 
functions," Z. Phys. B43 (1981) 119. 

[85] B. Baumann, "Noncanonical path and surface simulation," Nucl. Phys. B285 
(1987) 391. 

[86] B. A. Berg and T. Neuhaus, "Multicanonical algorithms for first order phase 
transitions," Phys. Lett. B267 (1991) 249. 

K. Kajantie and M. Laine and K. Rummukainen and M. Shaposhnikov, "The 
electroweak phase transition: A non-perturbative analysis," Nucl. Phys. B466 
(1996) 189, |hep-lat/9510020| . 



58 



[88] P. de Forcrand and O. Jahn, 11 SO (3) versus SU(2) lattice gauge theory," in 
Proceedings of the NATO workshop on "Confinement, Topology, and other 
Non-Perturbative Aspects of QCD", Stara Lesna, Feb. 2002, 2002. 
|arXiv : hep-lat/0205026[ 

[89] P. de Forcrand and O. Jahn, "Comparison of SO(3) and SU(2) lattice gauge 
theory," Nucl. Phys. B651 (2003) 125, lhep-lat/0211004| . 

[90] I. Montvay, "Embedding the M = 2 supersymmetric Yang-Mills theory in the 
adjoint Higgs- Yukawa model on the lattice," Phys. Lett. B344 (1995) 176, 
|arXiv : hep-lat/9410015 . 



[91] J. Wess and J. Bagger, Supersymmetry and supergravity. Princeton, 1992. 

[92] E. D'Hoker and D. Z. Freedman, "Supersymmetric gauge theories and the 
AdS/CFT correspondence," TASI 2001: Strings, Branes and Extra 
Dimensions, Boulder, Colorado, 3-29 Jun 2001, edited by S. S. Gubser and 



J. D. Lykken, River Edge, N.J., World Scientific (2004) |hep-th/0201253 



[93] M. F. L. Golterman and D. N. Petcher, "A local interactive lattice model with 
supersymmetry," Nucl. Phys. B319 (1989) 307. 

[94] J. Giedt, R. Koniuk, E. Poppitz and T. Yavin, "Less naive about 

supersymmetric lattice quantum mechanics," JEEP 0412 (2004) 033, 
|hep-lat/0410041| . 

[95] J. Giedt, "R-symmetry in the Q-exact (2,2) 2d lattice Wess-Zumino model," 
Nucl.Phys. B726 (2005) 210, lhep-lat/0507016| . 

[96] J. Giedt, "Symmetry and scaling in the Q-exact lattice (2,2) 2d Wess- Zumino 
model," PoS LAT2005 (2006) 270, |hep-lat/0508017| . 

[97] D. B. Kaplan, "A method for simulating chiral fermions on the lattice," Phys. 
Lett. B288 (1992) 342-347. 

[98] H. Neuberger, "Exactly massless quarks on the lattice," Phys. Lett. B417 
(1998) 141, |hep-lat/9707022| . 

[99] K. Fujikawa, M. Ishibashi, and H. Suzuki, "Ginsparg- Wilson operators and a 
no-go theorem," Phys. Lett. B538 (2002) 197, hep-lat/0202017 . 



59 



[100] M. Luscher, "Exact chiral symmetry on the lattice and the Ginsparg- Wilson 
relation," Phys. Lett. B428 (1998) 342, |hep-lat/98020l l . 

[101] Y. Kikukawa and H. Suzuki, "A local formulation of lattice Wess-Zumino 
model with exact U(1) R symmetry," JEEP 0502 (2005) 012, 
|hep r lat/0412042] . 

[102] W. Krauth, H. Nicolai, and M. Staudacher, "Monte Carlo approach to 
M-theory," Phys. Lett. B431 (1998) 31-41, |hep-th/9803iT7] . 

[103] T. Appelquist, G. T. Fleming, and E. T. Neil, "Lattice Study of the 

Conformal Window in QCD-like Theories," Phys. Rev. Lett. 100 (2008) 
171607, [ arXiv: 0712. 0609] . 

[104] Y. Shamir, B. Svetitsky, and T. DeGrand, "Zero of the discrete beta function 
in SU(3) lattice gauge theory with color sextet fermions," Phys. Rev. D78 
(2008) 031502, jarXiv: 0803 . 1707] . 

[105] T. Appelquist, G. T. Fleming, and E. T. Neil, "Lattice Study of Conformal 
Behavior in SU(3) Yang-Mills Theories," larXiv: 0901. 37661 

[106] J. Giedt, "Power-counting theorem for staggered fermions," Nucl. Phys. B782 
(2007) 134-158, lhep-lat/0606003j . 

[107] M. Luscher, R. Narayanan, P. Weisz, and U. Wolff, "The Schrodinger 

functional: A Renormalizable probe for nonAbelian gauge theories," Nucl. 
Phys. B384 (1992) 168-228, [hep-lat/9207009| . 

[108] M. Luscher, R. Sommer, P. Weisz, and U. Wolff, "A Precise determination of 
the running coupling in the SU(3) Yang-Mills theory," Nucl. Phys. B413 
(1994) 481-502, |hep-lat/9309005| . 

[109] B. Svetitsky, Y. Shamir, and T. DeGrand, "Nonperturbative infrared fixed 
point in sextet QCD," presented at 26th International Symposium on Lattice 
Field Theory (Lattice 2008), Williamsburg, Virginia, 14-20 Jul 2008. (2008) 
[arXTv: 0809. 2885] . 

[110] J. Greensite and B. Lautrup, "First order phase transition in four-dimensional 
SO(3) lattice gauge theory," Phys. Rev. Lett. 47 (1981) 9. 



60 



[Ill] G. Bhanot and M. Creutz, "Variant actions and phase structure in lattice 
gauge theory," Phys. Rev. D24 (1981) 3212. 

[112] S. Catterall, J. Giedt, F. Sannino, and J. Schneible, "Phase diagram of SU(2) 
with 2 flavors of dynamical adjoint quarks," JEEP 11 (2008) 009, 
[arXiv: 0807. 0792] . 

[113] T. DeGrand, Y. Shamir, and B. Svetitsky, "Phase structure of SU(3) gauge 
theory with two flavors of symmetric-representation fermions," Phys. Rev. 
D79 (2009) 034501, |arXiv:0812TT42 7]. 

[114] T. DeGrand, Y. Shamir, and B. Svetitsky, "Exploring the phase diagram of 
sextet QCD," larXiv: 0809 . 29531 

[115] G. D. Moore and K. Rummukainen, "Electroweak bubble nucleation, 
nonperturbatively," Phys. Rev. D63 (2001) 045002, |hep-ph/0009132| . 

[116] I. Horvath, A. D. Kennedy and S. Sint, "A new exact method for dynamical 
fermion computations with nonlocal action," Nucl. Phys. (Proc. Suppl.) B73 
(1999) 834. 

[117] M. A. Clark and A. D. Kennedy, "The RHMC algorithm for two flavors of 
dynamical staggered fermions," Nucl. Phys. (Proc. Suppl.) B129 (2004) 850. 

[118] M. A. Clark, A. D. Kennedy and Z. Sroczynski, "Exact 2+1 flavor RHMC 
simulations," Nucl. Phys. (Proc. Suppl.) 140 (2005) 835, 
[arXTv:hep-lat/0409133| . 

[119] D. Anselmi, D. Z. Freedman, M. T. Grisaru, and A. A. Johansen, 

"Universality of the operator product expansions of SCFT(4)," Phys. Lett. 
B394 (1997) 329-336, |hep -th/9608T25l . 

[120] F. Farchioni, C. Gebert, R. Kirchner, I. Montvay, A. Feo, G. Minister, 
T. Galla, A. Vladikas [DESY-Munster-Roma Collaboration], "The 
supersymmetric Ward identities on the lattice," Eur. Phys. J. C23 (2002) 
719, | arXiv : hep-lat/01 1 1008 . 

[121] J. Giedt, R. Brower, S. Catterall, G. T. Fleming, and P. Vranas, "Lattice 
super- Yang-Mills using domain wall fermions in the chiral limit," Phys. Rev. 
D79 (2009) 025015, |arXiv : 0810 . 5746 [hep-TatT] ■ 



61 



[122] J. Giedt, R. Brower, S. Catterall, G. T. Fleming, and P. Vranas, "Gluinos 
condensing at the CCNI: 4096 CPUs weigh in," in proceedings "Continuous 
Advances in QCD 2008," Fine Theoretical Physics Institute, University of 
Minnesota, Minneapolis, MN, May 15-18, 2008, ed. M. Peloso, World 
Scientific, Singapore [arXiv : 0807 . 2032 [hep-lat] ] . 

[123] J. Giedt, R. Brower, S. Catterall, G. T. Fleming, and P. Vranas, "Domain wall 
fermion lattice super- Yang-Mills," in proceedings "The XXVI International 
Symposium of Lattice Field Theory," College of William and Mary, 
Williamsburg, VA, July 14-19, 2008, PoS (LATTICE 2008) 053. 

[124] T. Blum et al., "Quenched lattice qcd with domain wall fermions and the 
chiral limit," Phys. Rev. D69 (2004) 074502, |arXiv: hep-lat /0007038] . 

[125] S. Ferrara and B. Zumino, "Supergauge invariant Yang-Mills theories," Nucl. 
Phys. B79 (1974) 413. 

[126] Y. Shamir, "Chiral fermions from lattice boundaries," Nucl. Phys. B406 
(1993) 90, |hep-lat/930300 5 . 

[127] F. Farchioni, G. Muenster and R. Peetz, "The volume source technique for 
flavor singlets: A second look," Eur. Phys. J. C38 (2004) 329, 
|hep r iat/0404004"] . 

[128] F. Farchioni and R. Peetz, "The low-lying mass spectrum of the N = 1 SU(2) 
SUSY Yang-Mills theory with Wilson fermions," Eur. Phys. J. C39 (2005) 
87, |hep-lat/0407036| . 

[129] "R. Peetz, "Spectrum of M = 1 Super Yang Mills Theory on the Lattice with 
a Light Gluino," doctoral dissertation, University of Miinster, Germany, 2003, 
available at: http : //deposit . ddb . de/cgi-bin/dokserv?idn=97018249X.." 

[130] "I. Montvay, Nucl. Phys. B83-84 (Proc. Suppl) (2000) 188; DESY-Miinster 
Collaboration (A. Feo et al.), Nucl. Phys. B83-84 (Proc. Suppl) (2000) 661; 
DESY-Miinster Collaboration (A. Feo et al), Nucl. Phys. B83-84 (Proc. 
Suppl) (2000) 670; DESY-Miinster Collaboration (R. Kirchner et al.), Phys. 
Lett. B446 (1999) 209; G. Koutsoumbas et. al. B63 (Proc. Suppl) (1998) 727; 
I. Montvay B63 (Proc. Suppl) (1998) 108; G. Koutsoumbas, I. Montvay, Phys. 



62 



Lett. B398 (1997) 130; I. Montvay, Nucl. Phys. B53 (Proc. Suppl) (1997) 
853; I. Montvay, Nucl. Phys. B466 (1996) 259.." 

[131] A. Donini, E. Gabrielli and M. B. Gavela, "Quenched super symmetry," Nucl. 
Phys. B546 (1999) 119, | arXiv : hep-th/9810127] . 

[132] A. Donini and M. Guagnelli, "Hybrid molecular dynamics for lattice 

supersymmetry," Phys. Lett. B383 (1996) 301, [|arXiv : hep-lat/96050l"0 . 

[133] Kamel Demmouche et al., "Spectrum of 4d M = 1 SYM on the lattice with 
light dynamical gluinos," in The XXVI International Symposium of Lattice 
Field Theory, College of William and Mary, Williamsburg, VA, July 14-19, 
2008, PoS (LATTICE 2008) 061. 

[134] M. Endres, "Numerical simulation of M = 1 supersymmetric Yang-Mills 
theory," in The XXVI International Symposium of Lattice Field Theory, 
College of William and Mary, Williamsburg, VA, July 14-19, 2008, PoS 
(LATTICE 2008) 025 [arXiv : 0810 . 0431] . 

[135] M. G. Endres, "Dynamical simulation of M = 1 supersymmetric Yang-Mills 
theory with domain wall fermions," [arXiv : 0902 . 4267 . 

[136] R. Sommer, "A new way to set the energy scale in lattice gauge theories and 
its applications to the static force and a s in SU(2) Yang-Mills theory," Nucl. 
Phys. B411 (1994) 839, |arXiv:hep-lat/9310022| . 

[137] RBC and UKQCD Collaboration, N. Christ, "Estimating domain wall 
fermion chiral symmetry breaking," PoS LAT2005 (2006) 345. 

[138] G. Martinelli, C. Pittori, C. T. Sachrajda, M. Testa and A. Vladikas, "A 
general method for nonperturbative renormalization of lattice operators," 
Nucl. Phys. B445 (1995) 81, |arXiv : hep-lat/941 lOTo] . 

[139] T. Blum et. al, "Non-perturbative renormalisation of domain wall fermions: 
Quark bilinears," Phys. Rev. D66 (2002) 014504, |arXiv:hep-lat/0102005l . 

[140] N. M. Davies, T. J. Hollowood, V. V. Khoze and M. P. Mattis, "Gluino 

condensate and magnetic monopoles in supersymmetric gluodynamics," Nucl. 
Phys. B559 (1999) 123, |arXiv : hep-th/9905015j . 



63 



[141] R. G. Edwards, U. M. Heller, and R. Narayanan, "Evidence for fractional 

topological charge in SU(2) pure Yang-Mills theory," Phys. Lett. B438 (1998) 
96-98, |hep-lat/9806011| . 

[142] J. Ambjorn, J. Giedt, and J. Greensite, "Vortex structure in 

Abelian-projected lattice gauge theory," Nucl. Phys. Proc. Suppl. 83 (2000) 
467-469, lhep-lat/9908017| . 

[143] J. Ambjorn, J. Giedt, and J. Greensite, "Vortex structure vs. monopole 
dominance in Abelian projected gauge theory," JEEP 02 (2000) 033, 
|hep-lat/990702i] . 

[144] L. Del Debbio, M. Faber, J. Giedt, J. Greensite, and S. Olejnik, "Detection of 
center vortices in the lattice Yang-Mills vacuum," Phys. Rev. D58 (1998) 
094501, hep-lat/9801027| . 

[145] J. Giedt and J. Greensite, "The abelianicity of cooled SU(2) lattice 



configurations," Phys. Rev. D55 (1997) 4484-4487, |hep-lat/9611002 



[146] R M. Vranas, "Gap domain wall fermions," Phys. Rev. D74 (2006) 034512, 



|arXiv:hep-lat /0606014 



[147] R. C. Brower, H. Neff and K. Orginos,, "Moebius fermions: Improved domain 
wall chiral fermions," Nucl. Phys. Proc. Suppl. 140 (2005) 686, |hep-lat/0|. 



64 



